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1. INTRODUCTION 


1 . 1 Statement of the Problem 

Multibladed pump inducers are used to increase the suction specific 
speed of centrifugal impellers for liquid rocket engines. Cavitation in 
rocket pumps is in general handled by inducers, which are secondary 
pumps designed to pass the required flow under cavitating conditions and 
produce enough head rise to permit the primary pump to operate 
cavitation free. Their main characteristics are very high solidity 
(ratio of the chord length of the blade to the blade spacing) and very 
low flow coefficient (ratio of inlet axial velocity to blade tip speed). 
While their specific geometry is dictated by cavitation requirements, 
the flow inside the long and narrow passages between the blades is sub- 
jected to major effects due to turbulence and viscosity, thus making the 
prediction of the flow behaviour in these inducers very difficult. 

The Department of Aerospace Engineering of The Pennsylvania State 
University has conducted investigation and analysis of the flow in pump 
inducers under NASA sponsorship since November, 19&3- It has designed 
and tested in air a four-bladed inducer model of three feet diameter 
(ref. **) . The investigation carried out on a three-bladed inducer of 
the same blade profile and flow characteristics is reported in this 
thesis . 

The secondary motions in these inducers are not confined to thin 
regions at the blades but extend over the entire cross section of the 
flow. The development of three dimensional turbulent boundary layers on 
the blade surfaces combined with the interaction between the suction and 
pressure surface boundary layers makes the flow truly three dimensional. 


J 

The principal objective of this thesis is to obtain a quantitative 
measurement of these secondary motions. It is hoped that a quantitative 
understanding of these secondary flows can serve the establishment of a 
theoretical model for the eventual analysis of the turbomachinery flow 
dominated by secondary fluid motions caused by viscosity and turbulence. 
In addition, attempts are made to establish the three dimensional 
inviscid flow through such inducers. Hence the purpose of this thesis 
has been: 

1) to develop a measuring technique using hot-wire probes to 
obtain the blade to blade variation of the three components of the 
velocity and turbulence at. the three bladed inducer exit. 

2) to predict the three dimensional inviscid flow in the inducer 
solving numerically the exact equations of motion. 

3) to approximately predict the flow at inducer exit using radial 
equilibrium and an empirical expression for the viscous losses. 

4) to investigate the axial flow velocity perturbation due to 
radially varying blade blockage. 

5) to carry out a complete flow survey using stationary probe at 
the exit of the inducer and a rotating probe inside the blade passage of 
the three-bladed inducer. 

6) to visualize the flow using ammonia streak technique and a tuft 
grid to derive respectively the limiting streamline angles on the blade 
surfaces and the maximum radial inward flow at mid passage for the 
three-bladed inducer. 

The performance of the three-bladed inducer is discussed thoroughly 
and compared with the performance of the four-bladed inducer previously 
studied in this department. 
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1.2 Previous Related Work 


A survey of the existing literature on axial flow Inducers reveals 
very little theoretical analysis for the prediction of the three 
dimensional flow characteristics. Although this particular case is 
related to the three dimensional flow problems of turbomachinery in 
general, the presence of large secondary flows caused by three 
dimensional boundary layers makes it essential to include the viscous 
effects in the solution of the flow field. Because of the lack of 
accurate information concerning the real flow processes and the 
complexity of the viscid equations of motion governing this flow, very 
little progress has been made toward a meaningful solution of the flow 
field. For the above reasons, most of the design and analysis of the 
Inducer fluid flow is based on conventional two dimensional methods. 

On the assumption of a fully developed turbulent flow along with 
assumed velocity profiles in both radial and mai ns t ream di rect ions of 
the channel and an empirically determined skin friction loss coefficient, 
an approximate solution using the simplified radial equilibrium equation 
in an integrated form has been obtained for a four-bladed inducer in 
reference 4 . This analysis leads to a good prediction of the outlet 
absolute tangential velocity. However, due to the lack of accurate 
information concerning the radial velocity profiles inside the blade 
passages, the axial velocities are predicted only qualitatively. A 
first step toward the study of the three dimensional boundary layer in 
such a channel is carried out by A. Jabbari (reference 3), who investi- 
gated the growth of the boundary layer on a single helical blade of 
constant' thickness rotating at the same speed as the inducer. A logical 
extension of this investigation to study the boundary layer growth in a 
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rotating channel to be undertaken by this department should lead to some 
valuable Information for an exact .prediction of the flow. 

An attempt to predict the exit head rise and flow coefficient for 
an 80° helical inducer by taking Into account the losses has been made 
by J. C. Montgomery (reference 11). This author uses the simplified 
radial equilibrium equation together with an arbitrary expression for the 
losses to predict the exit flow characteristics. Although his analysis 
predicts the large axial velocity gradient in the radial direction at 
exit, the procedure of trying arbitrary loss distributions to obtain the 
desired solution for the exit flow coefficient lacks the physical support 
It does not provide any specific method of relating the loss distribution 
to a given geometrical and flow characteristic of the machine. 

Among the attempts made toward a general solution of the equations 
governing the flow, the main contribution is due to P. Cooper and H. 

Bosch (ref. 2). Using a numerical procedure these authors solved the 
equations of motion iteratively in three dimensions for an exact solution 
of the flow field. This analysis, however, has its limitations and it is 
proposed to treat it in more detail in a following section. 

On the experimental aspects of the flow, several authors have 
studied the overall inducer performance with different configurations and 

: t 

inlet angles and tested in various fluids such as water, liquid hydrogen, 
or nitrogen under a wide range of flow parameters (see references 9, 10, 
11, 12, 14, 15, 16). However, in most cases these studies are very 
specific and deal with cavitation performance and overall efficiency 

I * ' 

rather than with a general and basic investigation of the flow phenomena 
in inducers. The cavitating and non-cavi tating performance of a 8^°, 

81° and 78° inducer under different flow coefficients has been studied by 


it 


Acosta (1). Acosta observes a deterioration in the radial distribution 
of the axial velocity and head rise at exit at a flow coefficient 
* = 0.07 under non-cavi tating condition, with a backflow near the hub 
and a sharp positive gradient in head rise near the tip. These results 
are similar to those obtained for the three-bladed inducer reported in 
this thesis. Acosta attributes the large departure from design values 
based on the Simplified Radial Equilibrium Equation (SRE) to strong 
three dimensional and viscous effects. Acosta also notices a decrease 


«? 

in efficiency at large solidity, which is confirmed by the results 

obtained for the b and 3 bladed inducers tested in this department. 

This effect can be attributed to the influence of blade blockage on the 

flow characteristics and an increase in viscous and turbulent mixing 

losses due to decrease in channel width. 

R. F. Soltis and D. A. Anderson (15) were led to similar 

conclusions while investigating the non-cavi tat i ng performance of a 78° 

axial inducer under various flow coefficients. Using the simplified 

radial equilibrium equation which relates the static head rise gradient 

and the tangential velocity: _ 

V 1 

9h .. V 9 

3r g Q r 

1 

They derived the outlet axial velocity profile using experimental values 

! 

for the total pressure and outlet flow angles. However, at their 

i 

station of measurement, located about 1/5 diameter downstream of the 
rotor, the radial velocities are likely to be very small and the flow 
nearly axisymetric since the wake diffusion in such inducers is very 
rapid. 

This analysis serves to establish that the flow is axisymmetric at 
small axial distances downstream of trailing edge (which may correspond 
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to large distance along the streamline path since the streamlines follow 
tightly wound helical path), but cannot account for the existence of 
large head rise near the tip. Similar observations are made by P. J. 
Mullan (9), P. R. Meng and R. D. Moore (10), J. C. Montgomery (11), W. 

M. Osborn (12). R. F. Soltis, D. C. Urasek and M. J. Miller (16) carried 
an investigation in tandem bladed inducer operated in water. Their 
results show an appreciable increase in efficiency and head rise compared 
to conventional inducers. Control of the boundary layer growth on the 
rear part of the blades in a tandem bladed Inducer should Improve the 
flow appreciably by reducing the interaction effects between pressure 
and suction surface boundary layers. 

The main conclusions of the various studies described above are: 

V 

1) The overall head rise coefficient increases when the operating 
flow coefficient decreases, especially in the tip region. 

'i 

2) The total head rise coefficient increases when the solidity of 
the blades is decreased in the practical range of solidities used in the 
i nducers . 

3) The radial distribution of outlet velocity tends to deteriorate 
when the flow coefficient is decreased. At low flow coefficients and for 
most inducer configurations, there is a large positive radial gradient in 
exit axial velocity with a backflow near the hub. This backflow is prp- 
babiy linked to a large redistribution of the flow as it leaves the 
trailing edge, the radial velocity decaying from presumably large posi- 
tive values inside the blade passage to zero as the flow becomes 
axisymmetric downstream. However, the radial velocity , at exit of 
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inaucers nas never been measured prior to this study. An attempt is 
made in this thesis to measure the radial velocity using a three 
dimensional hot-wire probe set-up.. 

1 . 3 Background of the Present Investigation 

As stated in section 1.2, the Department of Aerospace Engineering 
at The Pennsylvania State University has conducted the investigation 
of the flow in a four-bladed pump inducer. The three-biaded inducer 
reported in this thesis is of the same design and is described in 
section 3.1.1. 

The most important observations and conclusions obtained from four 
biaded inducer results are, (ref. A): 

1. At or near design flow coefficient, no back flow was observed at 
the inlet of the inducer. A separated region of flow exists near the 
hub at the discharge side of the inducer. This can.be explained quali- 
tatively on the basis of simplified radial equi 1 ibrium equation, when 
applied to the actual rather than the design head distribution. How- 
ever the actual region of separation is much smaller than that derived 
by the condition of radial equilibrium. 

The measurements indicate that this back flow region originates on 
the stationary hub just after it has left the rotating hub. Furthermore, 
this extent of back flow was observed to grow continuously downstream. 

The extent of back flow increases considerably at the inlet of tip 
region and outlet hub region at flow coefficients lower than the design 
value. 

2. The flow visualization experiments carried out near the exposed 
leading part of the inducer indicate that the radial motions in the 
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main flow (not under the influence of blade skin friction) are smaller 
where as the radial motions near the blade surfaces are considerable. 

3. The test inducer, designed approximately for uniform head distri- 
bution over its discharge area (assuming ideal flow) actually produces' 
a nonuniform head. Near the tip the actual head of the absolute flow 
was found to be 2 - 3 times that at hub and mid radius. This non-uni- ‘ 
form head distribution can be explained qualitatively by real fluid 
effects and is in agreement with the observations of other investigators. 

The flow survey carried out at the trailing edge and farther down- 
stream indicate that the radial distribution of the axial velocity 
component changes considerably between the two stations. The stag- 
nation pressure at the tip is found to decrease continuously as the flow 
travels downstream, whereas it remains essential ly constant at mid 
radius and hub. ! 

4. Measurement of the relative flow inside the blade passages (obtained 
with use of a rotating probe and pressure transfer device) indicate that 
a major portion of the total flow losses along the stream path occurs 
near the leading edge. This is probably'due to the presence of laminar 
flow on the blades near the leading edge and the consequent large radial 
flow and associated mixing effects near the tip. The flow losses at the 
tip were found to be very much higher than’ those at other radial 
locations. 

5. A complete flow survey of the relative flow inside the blade passage 
at the trailing edge revealed the presence of a loss core -located 
slightly inward from the blade tip. The radial movement inside the 
blade boundary layer, when encountered by the annulus wall tends to 
deflect towards the mid passage and then ‘radially inward. The mixing 
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ettects due to these secondary flows are -respons i b le for large losses 
observed experimentally. The measurements indicate the presence of 
large loss regions near the mid passage extending all the way from mid 
radius to the tip region, of the blades. 

6. The blade boundary layers were found _to be quite thin near the hub 
and mid radius, being thicker on suction surface than the pressure 
surface. 

The expected radial motions within the blade passages have been 
qualitatively confirmed by flow visualization experiments and appeared 
to be quite strong all along the blade length. 

The wall shear stress estimated from the boundary layer profiles 
inside the blade passage appeared to be much higher than that of an 
equivalent stationary channel having the same relative flow. 

7. At the trailing edge, the relative velocities averaged circum- 
ferentially over the passage had a maximum value near the mid radius. 

For .a third of the blade height from the tip the relative velocities 
were considerably less than their design values. 

The angles of the relative flow downstream of the trailing edge 
were found to be nearly uniform in the circumferential direction,, their 
deviation from the design values being greatest near the tip. 

8. A new frictionloss coefficient applicable to inducers operating in 
the range of flow coefficients $ = 0.065 to 0.2 was defined and derived 
fromthe inducer data available in the open literature. This empirical 
friction coefficient is found to increase exponentially towards the 
blade tip. For the Penn State four bladed inducer, the radial variation 
of frictional losses estimated from this newly derived empirical loss 
coefficient agreed closely with the measured values. 
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9. A circumferentially averaged . radial equilibrium equation was used 

to predict the relative and absolute tangential .velocities. The 
analysis was based on suitable assumptions for the radial and mainflow 
velocity profiles (based on the existing three dimensional turbulent 
boundary layer data available). The agreement between the theory and < 
experiment was reasonably good. . 1 

10. The axial velocities predicted, using the continuity and axial 

momentum equations with assumed radial velocity profiles and the ' 

derived tangential velocity distribution agreed qualitatively with the 
values measured from a stationary probe. 

11. Measurements carried out at several stations downstream (Ref. 5) 
indicated that the axial t veloci ty profiles under went marked change 
as the flow proceeded downstream. The back flow region near the hub 
developed after the flow had left the rotating hub and grew continuously 
as the flow proceeded downstream. 

12. Measurement of the instantaneous three dimensional velocities 
carried out using X configuration hot wire indicated that the radial ■ 
velocities ait the exit were of the same order of magnitude as the axial 
velocities over most of the radial locations. 

13. The limiting streamline angle of the flow near the wall measured by 
ammonia streak technique revealed that these angles are generally higher 
than those measured on a single blade (reference 3). The maximum 
deviation of the streamlines from the cylindrical surface (outside the 
boundary layer region) obtained by means of a tuft grid mounted at the 
exit of the inducer showed Substantial radial inward flow from mid 
radius to tip. 
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14. The. blade .to. blade variation of the- relative velocity derived from 
rotating probe measurements indicated that in addition to conventional 
boundary layer type of profiles near the' blades, a wake type of profile 
existed from mid radius- to tip. A qualitative analysis carried out 
indicated that the radial inward flow at mid passage at these radii was 1 
the fcause for the wake type of profile. 

The flow measurements for the four bladed inducer have indicated 
that the channel had a fully developed flow at most of the axial 
locations. In reference: 4, i t was anticipated that a decrease in the 
solidity of the blades might improve the performance of- the inducer by 
reducing the interaction of the pressure and suction surface blade f 

boundary layers near the tip. This was accomplished by- removing one of 
the blades. of the four bladed inducer and equally spacing the others, 
thus ^.decreasing the solidity by twenty-five percent. 

This thesis presents the performance of the three bladed inducer. 

- - . ! ' i 

1.4 Methods and Means of Investigation 

The complexity of the flow in inducers makes meaningful predictions 
extremely difficult. A flow model such as proposed in reference 4 based 
on .the assumption of fully developed turbulent flow and radial equilib- 
rium lead^ to a good approximation of the tangential velocity at exit, 
but does notipredict the axial velocity measured. This is due to lack 
of accurate information regarding the radial velocity distribution at 
exit. This approximate, viscid analysis is applied to the case of the 
three bladed inducer in section 2.3- 

In order to understand the extent of three dimensional effects, it 
is essential to have a knowledge of the inviscid flow field of the three 

1 1 
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dimensional flow taused by radially varying blade blockage, flow turning, 
finite hub to tip ratio and other effects. This is accomplished by using 
an existing numerical solution of the flow field proposed by Cooper and 
Bosch. This three dimensional method employs a numerical procedure to 
solve the equations of motion expressed in finite difference form and 
requires the use of computers with large storage capacity and fast ' 
computation time. This approach is discussed in section 2.2. 

In order to study the flow perturbations due to radial gradient of 
blade blockage, which is large in the case of the Penn State inducer, 
an analysis carried out by Lewis and Horlock is used. The predicted 
perturbations in axial velocity are found to be of appreciable magnitude. 

Thd final solution of the entire flow field, including viscous 
effects, can be carried out only by introducing some approximations in 
the three dimensional turbulent equations of motion governing this flow. 
Experiment is the bnly source which helps to make these approximations. 
Hence the following experimental program is unuertaken. It has been 
stressed that a good knowledge of the entire flow field is necessary for 
the design and prediction of inducer performance. Because of the three 
dimensional character of the flow within the blade passages, accurate 
measuring techniques are scarce. In this thesis, a method using a 
configuration of three hot wire sensors is developed to get the blade to 
blade distribution of the three components of the velocity at exit of 
the three bladed inducer. This technique is extended to the measurement 
of the local turbulence intensity and the method thus provides a good 
means of measuring the extent of secondary flow or the deviation from 
the design flow (section 3*3). 
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Finally, conventional experimental techniques such as flow 
visualization and the use of both stationary and rotating pitot tubes 
are used for measurement of the flow inside the blade passage and at 
several stations downstream of the inducer (see sections 3-2, 3.**, 3-5). 

i * 
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2. THEORETICAL ANALYSIS 


An accurate assessment of all significant perturbations in inviscld 
flow is essential for the final prediction of the flow including viscous 
effects. An important contribution towards the three dimensional 

i 

inviscld solution is made by Cooper and Bosch (reference 2). This 
solution is extended for the Penn State inducer in section 2.2 to i . 

determine the three dimensional inviscid perturbations. Although 
blockage is included in the above numerical method of analysis, the 
solution reveals little about the perturbations in meridional flow 
caused by various sources. The finite blade thickness and its 
variation in radial direction causes significant flow perturbations. 

This effect is considered in section 2.1. 

^ In section 2.3, an earlier analysis based on empirically deter- 
mined skin friction coefficient is used to predict tangential velocities 
at the exit of the three bladed inducer. 


2 . 1 Influence of Blade Blockage on Exit Velocity 

The axial velocities measured at the exit of inducers are generally 
found to be substantially different f rom> the design values. These 
departures are caused by (a) radially varying blade thickness - the ' 

blade blockage effect, (b) the radial flows inside the blade boundary 

. • * 

layers. ■ . . 

In order to find out the contribution of the blade blockage effect 

on axial velocity deviations the following analysis is carried out. 

The blockage factor is defined as 


6 = 1 - 


tn 

2irr 

14 


(1) 



t being the thickness of the blade in the tangential direction, 
n being the number of blades, 

■ r being .the local radius. 

which gives for the three bladed inducer (Fig. (2)): • 

6 = .800 at the hub 
.980 at the. t'ip 

and for the four bladed inducer: 

t . . .730 at the hub 

. .970 near the tip 

Such a variation of blade thickness from hub to casing imposes a non- 
uniform blockage upon the meridional flow; this should give rise to 
significant perturbations in axial velocity. 

R. I. Lewis and J. H. Horlock (reference 6) have developed an 
analogous source actuator theory accounting for the effects of vary- 
ing blade thickness, and have pointed out the significant effects on 
the meridional flow perturbations due to blade blockage. Following! 
Lewis and Horlock's theory, perturbations due to blade thickness are 
governed by 


' • 'divgrad<J>=S (2) 

<}> being'the scalar perturbation potential due to source d i s t r i but ion <and 
S the source distribution representing the finite blade thickness in'a 
uniform stream. In the case of known geometry, where an approximate* 1 
expression for the blockage factor 6(z,r) can be found, the source 
distribution S(z,r) per unit volume can be derived from the equation 
given by Lewis and Horlock: . : 


(r s<2 ’ r) ( z> 


(3) 
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where V z is the local undisturbed axial velocity without blockage, 
herein taken as unity. This expression is derived from continuity 
considerations between the blades for the incremental change in axial 
velocity TT. Knowing the blockage factor 6 from the blade geometry, the 
source function $(z,r) can be determined. Using this, equation (2) can 
then be solved. Lewis and Horlock's expression for axial velocity 
perturbation outside the blade passage is given by; 

V = V + v' 
z z z 



-(±)X. (z - z ) 
S(z s ,r) e 1 b 


dz 

s 


( 4 ) 





-(±)A. (z - 

S(z s ,r)e 


z s> 


drdz s ) 


r 


h 


t 

where 

v is the local resultant velocity 
z , 

v^ is the velocity perturbation due to blockage 
r^ and r^ are the hub and casing radius 
z g is the axial location at point of integration 
Zj and are respectively the leading and trailing 
edge axial coordinates. 

The. sign convention is + for z>Z 2 and - for z<z^. 


16 


Within the blade row the velocity perturbation is given by. 


v; = f ( j S(z g ,r) e 


-X, (z - z ) 


d z 


X.(z - z ) 

S(z ,r) e d z } 

s s 


Z 2 Tt 


2 2 
r t " r h 


(5) 

-X.(z - z ) 

r S(z g ,r) (1 - e ) dr dz g 


2 1 r h 


r S(z s> r) (1 


X, (z - z ) 


e s )dr d z^} 


z, r. 


Xj is determined by the condition of vanishing radial velocity at hub 
and casing: 

J 1 (X i r h } Y 1 (X 1 r t } " J 1 (A ! r t } Y ! (X 1 r h ) = 0 (6) 

and being Bessel functions of first and second kind. With values 

of r, = 1 and r = 2, this leads to 
n t 

X, = 3-1976 

It should be noted however that this analysis takes into account 
the perturbations due to blade blockage only, without allowing for 
perturbation due to annulus area changes. The latter effect is taken 
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into account in the definition of the undisturbed velocity V which is 

z 

based on the local velocity without blades. 

Equations k and 5 are integrated numerically for the geometry given 
in Fig. ( 2 ) for both four and three bladed inducers, for which the 
blockage factors can be approximated by: 

6 = 7.1 r 3 (5 - z) 2 - z ( 4.5 - z) 

7.1 r 3 (5 - z) 2 

for four bladed inducer. 

and for three bladed inducer: 

6 = p3 (5 - z 2 ) - 2 (*. 5 .- 2) (7) 

9 . A r 3 (r - z) 2 

(The axial coordinate is undimens ional ized with respect to tip radius, 

2 = z/r t ) 

The results plotted in Fig. ( 2 ) show that the perturbation In 
axial velocity due to blade blockage reaches 13% of the local unperturbed 
axial velocity for the four bladed inducer and 10% for the three bladed, 
at hub near the trailing edge. It is apparent that immediately beyond 
the trailing edge, the minimum in axial velocity occurs near the hub; in 
fact the varying blade blockage from hub to tip imposes a strong upward 
radial velocity especially near the hub which is responsible for a 
redistribution of the flow directed toward the tip. Figure ( 2 ) also 
shows that the gradient in axial velocity across the trailing edge is 
very large in both cases. It is to be noted that the qualitative nature 

t 

of the axial veloci ty i.e. decrease near the hub and increase towards 
the tip is predicted here by taking into account only the blockage effect 
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on the meridional flow (Fig. 19b). However there is no quantitative 
agreement between measured and predicted distributions of axial velocity. 
It is apparent from this result that an accurate prediction should 

r 

include a complete three dimensional solution with viscous effects 
Incorporated in the analysis. One of the useful conclusions of this 
analysis is that the blade blockage effect does contribute towards the 
decrease in axial velocity near the hub and its effect should be taken 
into account in the analysis for predicting the three dimensional flow. 

2.2 Numerical Solution of the Flow Field 

The availability of modern computers with large storage capacity 
and fast computation time greatly enlarge the possibility of solving the 
complete equations of motion. One of early attempts in this direction 
was carried out by Cooper and Bosch (reference 2) for the case of three 
dimensional inviscid flow through inducers. The author used these 
computer programmes, with some modifications to account for trailing 
edge boundary conditions, to predict the three dimensional inviscid 
flow through three bladed inducer. 

The quasi-three dimensional solution described in the following 
paragraph is used as a first approximation for the input data of the 

F 

three dimensional program in section 2.2.2. Because of the low conver- 
gence rate for a satisfactory solution of the flow field from the exact 
program, the initial conditions have to be as accurate as possible. The 
quasi-three dimensional solution, which predicts the flow quantities 
from what amounts to combining the axisymmetric solution with a blade to 
blade solution, is a better initial approximation than the values based 
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on one dimensional analysis. This approach results in considerable 
saving in the computer time required for the exact program. c v 

( ' ' ! t : • : r? •'-•••. / 

2.2.1 Quasi Three Dimensional Method ;• 

One of the computer programs developed by Cooper and Bosch derives 
the blade to blade average quantities using axisymmetric equations. 

These quantities are then used in a blade to blade solution using an 
integrated form of the scalar momentum equation in the tangential 
direction. A two dimensional meridional solution resulting from the 
streamline momentum equation together with the simple normal equilibrium 
equation is employed to determine the average quantities for the' velocity 
and pressure within the blade passage, 


v ^ 

dp 8 cos y* 
g odn 1 r 

L Is the friction loss term along the streamline 

m is the direction tangent to the streamline 

n 1 is the direction normal to the streamline 

Y 1 is the angle between the streamline and r coordinate. 


( 8 ) 

( 9 ) 


These quantities are then used in a blade to blade superimposed 
solution using an integrated form of the scalar momentum equation in the 
0 di recti on, 


iE = £ v <L { r u ) 

90 g m dm ' v 0' 
o 

together with the continuity equation. 


( 10 ) 
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These equations are solved numerically using a residual technique. 
The computer , program was somewhat modified to take into account the 
variation in blade thickness at each location within the grid of points 
representing the blade. The geometry and leading parameters are shown 
in Fig. (3). The grid is made of 26 points In the axial direction and 
7 points between hub and casing.' The solution of this program is used 
as input data for the exact three dimensional approach described in 
section 2.2.2. The accuracy criterium as defined in reference 2, 

t < e (fi 2 r ) £ (11) 

- P t g Q 

x being the streamline unbalance (or residual), is satisfied In about 

•J>C • ■* * . 

200 seconds on IBM 360-67 computer using a Fortran IV G compiler, with 

a value of 0.001 for e . 

P 


2.2.2. Exact Method of Solution 

This method developed by Cooper and Bosch solves the momentum and 
continuity equations iteratively in three dimensions for a finite grid 
of points representing the channel between the blades. The non linear 
partial differential equations governing the flow in a rotating polar 
coordinate system r, 8, z are: 


r momentum 


g » 9W W Q 3W 9W . 

— + W -*£+-£ r-llu + 

p 9r r 9r r 30 z 3z r 0 


rfi) ‘ 


+ F = 0 
r 


0 momentum 

g. , n 9W fl w 9W 3W w w 
— !§ +W + — TiT+ w 3-9 + -£-9 + 2 w n + F fl = 0 

p 90 r 9r r 30 z 9z r r 0 
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z momentum 

g . 3W W* 3W 3W 

-°|£+w * + -i — i + W — i+F 

p T5z W r a r r 30 z 3z 2 


( 12 ) 


continuity i 

w aw , aw Q aw . _ w Q . . 

_n + c. + J. §. +' —2. + J. ( W M + J.M+i w i£) _ o 

r 38 r 30 3z p ' r 3r r 30 z 3z' 


These equations, where is the angular velocity of the inducer, 

W , W_, W the relative velocities in radial, tangential, and axial 
r a z 

directions respectively, are expressed in a finite difference form and 

t 

rearranged to give four residuals which are reduced to zero by a 
relaxation procedure. A complete description of the numerical technique 
and the computer program can be found in reference 2. 

This method of solution has been applied by the author with some 

V i 

modifications. The flow is assumed to be incompressible throughout thus 
eliminating the last three terms Inside the bracket in the continuity 
equation above. A grid of 7 ( X 7 x 26 is chosen to represent the blade 
passage, and the initial values at each of those grid points are derived 
from the approximate solution results, described in section 2.2.1. As 
shown in Fig. (3), 7 stations are chosen in both the radial and 

l > 

tangential directions. The boundary condition to be satisfied on the 
hub, annulus walls, and the blade surfaces is 

V? • n = 0 ( 13 ) 

where is the direction normal to the channel boundaries. 

Among the 26 equally spaced axial stations, the first station 
corresponds to the upstream through flow boundary where the initial 
conditions are applied. For the boundary value problem to be consistent, 
these initial upstream conditions must specify the three components of 
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the velocity and the pressure, and the tangential velocity on the 
second axial station which thus defines the whirl at inlet 'of the 

• ■ i 

inducer, herein taken as zero. The four last stations correspond to the 
downstream through flow boundary, and extends about one fifth of the 
chord length downstream of the trailing edge. With the condition (13) 
to be satisfied on these stagnation stream surfaces, the set of boundary 
conditions for this problem is complete. The program has been suitably 
modified for use in IBM 360-67 system at The Pennsylvania State 
University. The solution obtained from this program was checked against 
a known solution for the case of a paddle wheel channel similar to that 
carried out by Cooper and Bosch (reference 2). 

a) Fi rst Run 

The program was first run for twenty-four cycles for the Penn State 
inducer, with the geometry shown in Fig. (2). The design air inlet 
angles at the leading edge are shown on Fig. (6). Normally, as sug- 
gested by P.' Cooper and H. Bosch (2), the upstream through flow boundary 
should be applied about one channel width upstream of the leading edge 

t . 

to take into account the upstream effects of the channel. Since the 
Penn State inducer is designed for light loading on the inlet portion of 
the blades, the upstream perturbation due to leading edge loading can be 
considered negligible. 

Furthermore, four axial stations have been chosen to represent the 
downstream extension of the blades, accounting for 1/5 of the blade 
chord length. 

The extension of the stagnation stream surfaces downstream with the 
condition that no flow crosses them' sat i sfies the requirement that the 
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flow be uniformly periodic at exit with spacing of 2i\/n, n being the 
number of blades. The extension of these surfaces, together with the 
manual method of, adjusting them to unload them,. has been discussed by 
P. Cooper and H. Bosch. 

The blade static pressures obtained by this program are shown on 
Fig. (5a, 5b, 5c) at hub, mid radius, and tip respectively. The 
negative pressures and gradient of pressure observed near the leading 
edge ar;e consistent with hub and wall static pressure probe measurement. 
In view of the results shown in Fig. (5a, 5b, 5c) > it can be said that 
the condition of closure of the blade pressure diagram at the trailing 
edge is far from being satisfied. As indicated by Cooper and Bosch, 
probably more stations would be required downstream for the stagnation 
stream surfaces to satisfy the closure condition. This method requires 
considerable amount of computer time, and hence is not the most efficient 
method. Cooper and Bosch also tried to satisfy this condition by forcing 
the pressures according to the equation, 

P = £_ ( r 2 n 2 + w ] 2 - w p 2 w 0 2 - w z 2 ) (1 k) 

9 o 

Wj being the relative velocity at inlet at radius r. This method is 
subject to criticism since equation (1**) is valid only alongj a stream- 
line and not at any radius. 

The main effort of the work reported here has been to try more 
efficient methods to satisfy the Kutta- Joukowsk i condition at the 
trai 1 ing edge. On Dr. Wislicenus 1 suggestion, a first try with forced 
stagnation pressure at the trailing edge (relative velocity equal to 
zero) has been attempted without success. The coarseness of the grid 
along with input problems did not allow for a detailed specification of 


2k 


I 



1 


the trailing edge. According to the sketch below, by locating and 
fixing the stagnation point on the trailing edge, it is possible to 
simulate real fluid effects. However, the grid of points has to 
represent in detail the trailing edge, and consequently has to be very 
dense in order for the large velocity gradients induced by this method 



A second attempt was made by setting the exit direction of the flow 
free, that is the boundary condition 

W • n = 0 (15) 

is no longer required on the downstream stagnation stream surfaces. For 
the mathematical problem to be consistent, this condition is replaced by 
setting the pressure on these downstream extensions of the blades equal 
to the averagd between pressure and suction surfaces at that particular 
radius : 


P 


p + p 

Li j e. 


( 16 ) 


Mathematically, this condition is equivalent to the boundary condition 
(15). Hence, according to the discussion given by P. Cooper and H. 
Bosch (2), the number of independent discrete variables is the same as 
in the previous case, and the solution to this problem exists and is 
unique. With this modification the condition of closing the blade 
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pressure diagram is perfectly fulfilled, Fig . (5a, 5b, 5c). The exit ( 
flow direction is then determined byi both axial and tangential 3 

velocities, which satisfy the equations of motion and continuity. 

Although the frequent manual manipulation of the downstream 
extensions of the blades is eliminated by this method, it should be - 
stressed • that the new condition specified by equation (16) is valid only • 
on the downstream stagnation surfaces. Because the static pressure 
coefficient is not necessarily constant i;n the tangential direction at 
any fixed radius, this condition would not prevent a tangenti-al pressure 
gradient to exist across the stagnation stream surfaces. Consequently, 
these surfaces have to be set up at the beginning of the numerical 
process from the exi t idevi at i on angles predicted by the design or any 
approximate method. As the computation proceeds, they can be readjusted 
after a few cycles from the exi t deviation angles derived'from>the 
program. With such a procedure, a method for automat i cal ly adj us t i ng 
the stagnation stream surfaces could be relatively eas i ly i ncorporated 
in the program. 

Although more cycles would have been required for the sol ut ion.vto ; "... 1 
stabilize, the results of Fig. (5a, 5b, 5c) showlthat the Kutta- 
Joukowski condition for the blade pressure distribution is satisfied. 

Also plotted on Fig. (5a, 5b, 5c) is the experimental rise in static 
pressure through the inducer. This experimental data is obtained 
through static ta'ps on the casing and holes drill.ed on the hub midway t 
between the blades and from static pressure taps on the blade surface. s - 
The agreement between predicted static pressure and experiment i,s good, r 
at all radii. The average flow properties obtained by this method-at 
the inducer exit are shown plotted in Fig. (A). The radial ; di stribut ion -■ 
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of axial velocity shows a trend similar to that of the experiment 
(stationary probe measurements). The discrepancy in magnitude between 
the predicted and experimental results is due to radial velocities 
induced by viscous effects near the blade surfaces. 

The absolute tangential velocity, although substantially higher 
than the design value, is approximately; constant from hub to mid radius 
and then decreases near the tip. Wislicenus means treaml Ine method 
(reference 1 7)» on which the design of the Penn State inducer, is based, 
predicts deviation angles of the order of 18° at the root and 5° at the 
tip. These large deviations are probably due to the empirical equation 
used by Wislicenus for deriving the deviation of the camber line from 
the mean streamline. These empirically determined deviations are based 
on NASA cascade tests carried out at moderate solidities. This method 
is yet to be perfected to make it valid for all types of turbomachinery. 
Hence the discrepancy between the results of the three dimensional 
program and the design values may be due to; 

1) The three dimensional effects such as blade blockage and the 
large^ci rcumferential extant of the blades* 

2) The inaccuracy of the means treaml ine method, 

3) The crudeness of the method adopted in this thesis to satisfy 

the Kutta-Joukowski condition. > 

The outlet deviation angles are slightly lowerthan those of the 
design as shown in Fig. (6). Since the stagger angles of the blades are 
large, a slight change in deviation angles results in considerable 
Change In absolute tangential velocities. This accounts for the large 
discrepancy between the predicted and design values of the absolute 
tangential veloci ty* shown in Fig. (4).' 
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The trend observed for the radial velocity distribution is 
probably due to blade blockage effect. The experimental' values of the l 
radial velocity derived from hot wire measurements, Fig. (12a), are very 
much higher than the predicted ones. Here also, the importance of the 
viscous effects which are neglected in the present analysis probably 
accounts for much of the discrepancy observed. The radial distribution 
of the static pressure coefficient shown in Fig. (k) Is consistent with 
the absolute tangential velocities and is higher than the design which 
assumes larger flow deviation angles at the trailing edge. 

In summary, it can be concluded that the program with the above 
modifications for pressure distribution provides a good means for 
predicting the flow in axial inducers. The accuracy depends on the 
number of grid points which can reasonably be chosen. An ultimate step 
would be to include real effects to the solution through ,the loss terms 
in equation (12). A loss subroutine was constructed for this analysis 
but has not been included in any additional run because of cost problems. 

The program was used with a Fortran IV H compiler using an 
optimization procedure which shortens the time required for repetitive 
calculations. It was executed on an IBM 360-67 approximately ten times 
faster than the Univac 1107 system used by P. Cooper and H. Bosch. 

2 . 3 Approximate Viscid Solution 

The two inviscid analyses carried out so far do not predict 
accurately the flow in axial inducers. It is apparent that the viscosity 
of the fluid in the long and narrow passages between the blades accounts 
for major departure of the flow characteristics from the design values, 
especially near the tip. In addition, the secondary flows induced by 
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the interaction between the pressure and suction surface boundary layers 
near the tip produce complex flows and Increase the flow losses. 

Although the prediction of these radial velocities and the flow losses 
Is extremely difficult, a first attempt to Include viscous effects has 
been made by B. , Lakshmi narayana (reference 4). In view of the great 
difficulties Involved In solving the complete equations of motion of the 
turbulent flow inside the passage, B. Lakshmi narayana sought to account 
for the viscous effects In an approximate way by using an empirical pres 
sure loss coefficient derived from various tests carried out at the NASA 
Lewis Research Center, M.l.T. Gas Turbine Lab., and TRW Cleveland 
(references 18, 19, 20 , 9 ). This analysis takes Into account the three 
dimensional nature of the flow and energy losses due to friction in 
solving continuity, axial momentum and circumferentially averaged radial 
equilibrium equations. The Inducer passages are assumed to have fully 
developed turbulent flow. The assumption for the velocity profiles In 
both tangential and radial directions are based on the present state of 
the knowledge of the three dimensional boundary layer. Fig. (7). 

An empirical friction loss coefficient applicable to inducers 
operating at low flow coefficients is derived from tests reported in 
references (18, • 19, 20 , 9 )• The friction loss coefficient in the 
following equation is plotted in Fig. (7). 


^ 9 o^Loss _ , ^ht 1 \ 

-75— — twi; V 

t n 


where , Fig . (7) : 


H. is the head loss due to friction 
Loss 


ht 


is the ratio of hub and tip radii 
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R is the Reynolds number Wd,/v 
n n 

Att r 

d, is the hydraulic mean diameter » — — cos B 
n n 

L q is the length of the channel 

a is the limiting streamline angle with m coordinate 

3 is the angle of the meridional component of the flow with the 

inducer axis 


The radial equilibrium equation, which incorporates all these assumptions 
and is valid for any noncavi tating inducer, is derived from equation (12). 
The resulting radial equilibrium equation is, 


n 97 c _ 9W fi _ o 9ot _ 9W e 

• -47 w 8 (tan “ 3^- + “e sec a 7T ) * '• 0IS w e ' S7~ 

cos S 


W* 2 X D R. . L ..2 

- 1-015 -4 (-U frr>-» 

R n h 

where x is the coordinate parallel to the blade camber line and lying 
on a cylindrical surface. Wq and W are defined on Fig. (7). 


The basic assumptions made in deriving this equation are: 


1 . The loss term F r in the radial momentum equation Is small. 

2. The radial gradient in axial velocity is very small at any 


particular radius. 

3. Tha flow properties are symmetrical about the mid passage 
(see Fig. (7) ) . 

i*. The relative velocity W at any radius is parallel to the mean 
blade camber line at that radius. 
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5. Wg varies linearly with x, the distance along the blade, from 
leading edge to trailing edge. 

6. The angle a between the limiting streamline and W varies 

m . . 

linearly from zero at the leading edge to the maximum' measured 
value at the trailing edge. 

In the case of the Penn State inducer, equation (17) reduces to: 

2 

3W W 

V 3 “ 0-015 + 2 F (r) ) + (1.015 + rF 1 (r) - £ (r)) 


- 2fl W Q (1 - £ (r)) 

wi th 

C(r) * - ■ (tan a + a sec 2 a) 
si n B 


(18) 


^( r ) o 0-03 tana 
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F(r) * j 


X R 
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ht 
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A R is the friction loss factor, derived in reference 4 and shown in 
Fig. (7)- Its validity has been proved for the range of flow coef- 
ficients ♦ = 0.065 to 0.2. F ( r) is shown plotted on Fig. (7) for the 
three b laded inducer. 

Equation (18) is solved for the three bladed inducer with the loss 
coefficient A^ given by B. lakshminarayana. The values of the functions 
C(r), C ( r) are the same as for the four bladed Inducer, reference 4. 

The boundary condition assumed for the solution of the differential 
equation (18) is that the tangential velocity Wg at hub is the same as 
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that of design. The predicted value of the absolute tangential velocity 
Vg is shown plotted versus radius in Fig. (8). It can be seen that the 
agreement between predicted. values and values measured from a. stationary 
pitot tube at station k, Fig. (10), is good. The steep rise in 
tangential velocity near the tip is confirmed by this analysis. 

It is quite clear from. this analysis that ifthe precise nature of 
velocity profiles and shear stresses is known, it is possible to 
predict the flow accurately. Theoretical and experimental investigation 
of the boundary layer in a rotating helical channel, currently under 
investigation in this department, should provide valuable information 
for the accurate prediction of the flow properties in such low flow 
coefficient, high solidity inducers. 
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3.' .EXPERIMENTAL PROGRAM AND RESULTS 


In carrying but the following experiments, the main goal is to 
investigate as completely as possible the flow characteristics such as 
time averaged and instantaneous veloci t ies , s tat i c and total heads of 
the relative and absolute flow. The importance of this experimental 
data for a better- understandi ng and prediction of the flow has been 
pointed out in section 1.1. It is to be stressed here that this study 
deals only with the noncavi tati ng performance of the inducer, aimed at 
understanding the basic flow phenomena. 

3.1. Apparatus Used in the Experiment 

3.1.1 Three B laded Inducer 

The experimental investigation is carried out in a three foot 
diameter model designed and built in the laboratory. The Inducer test 
facility is shown in Fig.(l). This model was originally studied with A 
blades, and later with three equally spaced blades. The blades were 
designed by G. F. Wislicenus and B. Lakshminarayana following the mean 
streamline method (reference 17). The design blade and flow angles are 
specified on Fig. (6). The inducer model has a rotational speed of A50 
RPM, giving a Reynolds number of 7 x 10^ based on tip radius. The RPM 
is determined to an accuracy of one tenth of RPM by means of a photocell, 
disk and electronic counter. Important parameters of the inducer are 
as follows: 
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Number of blades 


3 


Hub/Tip ratio (outlet) 3 0.5 

Tip diameter = 36. 0 inches 

Design value of flow 

coefficient (inlet) = 0.065 

Blade tip velocity = 71 feet/second 

Blade chord at 

Tip section = 82.96 inches 

Mid-span = 63*18 inches 

Hub = 49.94 inches 

Holes on the blade suction and pressure surfaces are connected to 
the metallic tubes embedded in the blade and are used for the limiting 
streamline angle measurements as well as the blade static pressure 
measurements. Figure (10) shows the location of the stations of measure- 
ment along the inducer axis. Stations 1 and 2 are located Inside the 
blade passage and are used for rotating probe measurements. Station 2a 
is located slightly outside the blade passage, distant of .25 inches 
from the blade trailing edge. Hot wire measurements are taken at station 
2a. 

3.1.2. Conventional Probes 

The total pressure probe used for stationary and rotating probe 
measurements is described in reference 8. The outside and inside 
diameters of the probe are respectively 0.125 and 0.055 inches. This 
combination gives a good sensitivity for a Mach number up to 0.045 
which is within the range of velocities encountered in the present 
investigation. The static pressure probe used in the stationary probe 
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measurements is standard, with four static holes each having a diameter 
of 0.010 inches and located 0.20 inches from the probe tip. 

For the measurement of the flow angles at exit a tuft probe made 
of teflon material is used. The same type of tuft is used for the 
qualitative measurement of the radial inward flow at the exit. 

3 . 1 . 3 . Appara tus for Instantaneous Velocity Measurements, Fig, (ila) 

Hot Wire Anemometer and Probes : 

Two dual channel hot wire constant-temperature anemometers designed 
and developed by the turbulence group at The Pennsylvania State 
University are used for this experimentation. Their original circuitry 
is shown on Fig. (lib), which is slightly modified for higher output 
voltage level. The Thermo Systems sensors used are made of tungsten hot 
wires of 3 y diameter and a length to diameter ratio of 700 . The probe 
configuration and mounting mechanism are shown on Fig. ( 9 ). 

An x configuration probe whose wires are aligned in the tangential 
and axial directions and a single wire probe aligned In the radial 
direction are used for these measurements. They are mounted on a 
traverse gear such that position error is negligibly small (Fig. (9)). 

Peripheric Instrumentation 

A very loW turbulence calibration tunnel built in this laboratory 
is used for calibration of the hot wires. This wind tunnel has a cross 
section of 2 x 2 inches and operates within the range of air velocities 
of 0 to 300 feet per second. 

A precision mi cromanometer manufactured by Flow Corporation is used 
to obtain the pressure measurements in the investigation. This manometer 
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has a minimum reading scale of 0.0001 Inches. Its accuracy is ±0.09 
per cent of the minimum static pressure measured in this particular 
experimental program. 

Digital and RMS voltmeters manufactured by Hewlett Packard and 
Ballantine Laboratory Inc. respectively are used for instantaneous 
velocity measurements. The data is recorded on a CEC data tape analog 
tape recorder at a speed of 3-75 i ps and displayed at 30 ips in an 
analog to digital converter (ADC) unit (Pastoriza Electronics, Inc.). 

The signals are vi sual 1 zed on a 4 channel storage Tektronix 600 
oscilloscope. Finally, the digital tapes are processed by the 

> 

University IBM 360-67 and IBM 1401 (ca 1 comp lotter) computers. 

A block diagram of the instrumentation used for the Instantaneous 
velocity measurements is shown in Fig.' ( 1 1 a ) . 

3.1.4. Pressure Transfer Device (PTD) 

A one channel pressure transfer device Is used for the rotating 
probe measurements and is described in details in reference 3. This 
device is mounted inside the hub nose cone of the model inducer and 
connected to a micromanometer. A similar ammonia transfer device is 
used to transfer ammonia from a static container to the location at 

which limiting streamline angles are to be measured. 

i i 

3.2 Flow Survey at the Exit Using Stationary Probes 

A complete flow survey at the exit of the inducer is carried out 
at stations 3, 4, 5 and 6 shown in Fig. (10). Total pressure and angles 
of the absolute flow are measured using total and tuft probes. The 
static pressures are measured at the hub and annulus walls at these 
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stations.- The rad i al::.d i s tr i bu t ion of ' the static head coefficient is 
derived by„) using the simplified radial equilibrium equation near the 
walls, 

V ^ 

•i lE = ’JL ■ 

p 3r r 

* it . 

This distribution has been compared with the one derived from a static 
probe and both were found to be consistent. These results are discussed 
and compared with similar measurements made for the four bladed inducer 

previously studied in this department and reported in references 4 and 5- 

1 ■ t 

: » * 

Stagnation and Static Pressure Coefficients 

The radial distribution of stagnation pressure coefficient (ip ) 

t ! 

and static pressure coefficient (i at various stations are shown 
plotted and compared with the higher solidity four bladed inducer 
resul ts in Fig . (19a). 

; ‘ ; i 

The static pressure coefficients are found to be consistently 
higher than those of the four bladed inducer (about 35 percent) at all 
radii. Surprisingly the radial gradient of static pressures remains the 
same in both cases. A slight decrease in the gradient is observed as 
the flow proceeds downstream. Fig. (19a). 

The stagnation pressure coefficients are found to be about 15 to 
20 percent higher than those of the four bladed inducer except near the 
tip. Furthermore, it can be observed by comparing distribution at 
station (3) and, (4) that the mixing losses downstream of the present 
inducer still exists but in fact are appreciably lower than those of the 
four bladed inducer. Even though the steep increase in the values of 
<p towards the tip is there even in the present case, this steep rise 
is confined to a smaller portion of the blade height near the tip than 
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tne one measured in the tour bladed inducer. There is appreciable 
improvement in the stagnation head distribution from hub to mid radius. 
No appreciable improvement in the head distribution is observed near the 
tip, and hence additional modifications in the form of boundary layer 
control or redesign are necessary to obtain an acceptable energy 
distribution. 

Radial Distribution of Axial Velocity 

The radial distribution of axial velocities at station 3, **, 5, 6 
is shown plotted in Figure (19b). 

The back flow is present at all downstream stations except at 
station 3- At station 3 the radial distribution of axial velocity is 
more uniform than that observed in the' four bladed Inducer and has 
appreciably higher value at mid radius. This indicates that the blade 
blockage effect which is higher in the case of a four bladed inducer has 
appreciable effect on the axial velocity distribution at the trailing 
edge of the blades. The axial velocity profiles show a marked 
deterioration from that observed in the four bladed inducer at all other 
stations. The extent of back flow region near the hub is appreciably 
higher than that of the four bladed inducer. 

t 

The predicted axial velocities from the exact analysis described 
in section 2.2.2 are shown plotted and compared with experimental values 
for the three bladed inducer in Fig. (4) . This plot shows the influence 
of the blade blockage. It is clear, from the observed discrepancy 
between the predicted (inviscid) and measured values, that it is 
essential to include the radial velocity induced by viscous effects for 
accurate prediction of axial velocity. 
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Radial Distribution of Absolute Tangential Velocities 
The radial variation of absolute tangential velocities is plotted 
in Figure (19c) . 

At stations 4 and 5 the tangential velocities from hub to mid radius 
are higher than those observed in the four bladed inducer whereas from 
mid radius to tip a slight decrease is observed. To some extent there 
Is an Improvement in the Vg distribution. These measurements are 
consistent with stagnation head rise coefficients described earlier. 

Even though the average increase in tangenttal velocities is small, there 
Is considerable (about 15% to 30%) Increase in stagnation head coef- 
ficients at these locations. The ratio of the actual head to the Euler 
head is a measure of the hydrodynamic efficiency. It is clear that 
there is appreciable decrease in the frictional effects when the 


solidity is decreased. The hydrodynamic efficiencies are generally 
higher in the case of three bladed inducer. Further improvement in the 
performance is possible only by boundary layer control. 


The measured is shown compared with values predicted from 
scid theory and approximate viscid theory on figures (4) and (8) 


respective ly 


Hub and Annulus Wall Static Pressures 

Wall static pressure coefficients measured on the hub wall with the 
aid of a pressure transfer device and on the annulus wall are plotted in 
Fig. (20) and (21) respectively. 

The values of pressure coefficient measured at all axial stations 
are considerably higher . .than those of a four bladed inducer. Further- 
more, the negative radial . pressure gradient observed at some of the axial 
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stations of a four bladed inducer has almost disappeared. The magni- 
tude of the negative pressures observed near the leading edge is , 

considerably reduced in the present case. .These results again confirm 
substantial improvement in the flow behavior when the solidity is 
decreased. • 

3.3. Instantaneous Velocity Measurements at Inducer Exit 

3 . 3.1- Abstract of the Method , 

In order to obtains better understanding of the flow phenomena- 
in axial flow inducers, an attempt is made to measure all thr;ee 1 
components of the three dimensional instantaneous velocity, at the ; 
inducer exit using a three hot-wire probe and the -associate equipment 
described in section 3. 1 - 3 » Fig. (11a). As pointed out earlier in 1 
Chapter 2, the radial velocity components are the most difficult to 
predict and measure. However, from a good knowledge of their magnitude 
and distribution, considerable improvement can be made in the prediction 
of the flow. 

The block diagram describing the method of measurement and The 
equipment required is shown on Fig. (11a). Those measurements are 'taken 
at station 2a, Fig. (10), which is located a quarter of an inch from the 
trailing edge of the blades. Since the radial velocity decays rapidly • 
downstream of the trailing edge, it is essential to 'carry out these 1 
measurements as close to the trailing edge as practicable. Since, the ; 
three ' components of the velocity constitute three unknowns, 'it is '• 

necessary to use three sensing wires to derive the magnitudes of the 
three velocity components (Vg or Wg, V r> V . .The probe traverse 


mechanism set up on the annulus wall is fixed tangentially, but is 
movable in the axial and radial directions (Fig. (9)). 

The three signals picked up from the hot wire anemometers are 
memorized in a magnetic tape, then digitized on an analog 1 to digital 
converter unit, (ADC), to be finally processed on an IBM 360-67 compute 
Fig. (11a). 

Considering the 1 fixed frame of reference r, 0, z, each one of the 
three wires is aligned with one of the coordinate axes. The three non- 
linearized voltage signals picked up by the probe are called Ej, 

Since the angle between the wire in 0 direction and the resultant 
velocity is very small, a correction to the deviation of the cosine law 
is used. According to King's law, the signals picked .up by the three 
sensors are, 


E 1 2 = A I + b i <v ' /2 

E 2 2 - A 2 + B 2 <V' /2 

E 3 2=A 3 +B 3 (V r6 ),/2 


(19) 


The sens i 1 1 vi ty .of the hot wi re to the flow direction is generally 

assumed to be given by the normal component or cosine law, whereby the 

effective cooling velocity is obtained from the magnitude of the vector 

component of the free stream normal to the axis of the wire. However, 

it has been proved, that i the f i ni te values of the length to diameter 

ratio,- 1/d, of the hot wires produce deviations from the cosine law. 

According to the sketch below, the effective cooling velocity in a 

free s.tream V is ■ , 
o’* 


a being the angle of incidence between the free stream direction and the 
normal to the axis of the wire, and k is an experimentally determined' 
parameter reflecting deviations from the cosine law. This deviation is 
important in the case of the wire sensing the resultant of the radial 


and axial velocity as shown in the following sketch: 



With this correction the velocity sensed by each wire is given by: 
V Z0 ■ ( v e 2 * " 2 2 >’ /2 
V <V 2 2 *Vr 2 + k 2 V e 2 > l/2 

v r e • (V r 2 ♦ v a 2 > ,/2 


leading to the following equations: 


V 2 = (_£i 
0 v 2 


1 -k 2 /2 
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W. H. Schwarz and C. A. Friehe (13) have carried measurements for the 

deviation from the cosine law which give for this case a correction 

factor of k = 0.26 for 1/d = 700. With this correction, the coefficients 

in the above equations are: 
k 2 

1 - = .966 


1 - k = .932 


Since the radial and axial velocities are an order of magnitude 
lower than the tangential velocity these corrections are not small.- 
The relative velocities are then given by, 



Wg = Vq - Qr (21) 

W = V 
z z 

The signal processing is carried out In five steps: 

1) The signals from the tape recorder are digitized in the ADC unit 
as input data for the computer program reproduced in Appendix 2. 

2) The digitized values of the signals are arranged in two 
dimensional arrays, this being dictated by the averaging method discussed 
in the following section. 

3) All digitized values are then linearized individually according 
to equation (19), with the parameters Aj and Bj accurately determined by 
probe calibration. 

4) Then the values of the resultant velocities are obtained from 
equations (20) and (21) to give the three components of the instantaneous 
velocity, and averaged over consecutive blade passages (see section 
3.3.2). 

5) Finally the turbulent intensity is computed for each component 
(see section 3*3.3). 

A computer program wr i tten for processing the digitized data is 
given in Appendix 2. 
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3 . 3 . 2 . Averaging Procedure 


The three signals are functions of time with a periodicity equal to 
three times the rotative speed of the machine. The only requirement for 
the condition that each sample. of the signal corresponds to the same 0 
location within the blade passage is that the sampling frequency be a 
multiple of the periodicity of the signal. .In other words, the following 
relationship has to be satisfied between the signal frequency (3ft), the 
number of samples n^ within the blade passage and the sampling frequency 
f of the ADC unit: 

n x 3ft = f - - (22) 

o 

The ADC unit used in the Department of Aerospace Engineering 
of The Pennsylvania State University has a fixed digitizing speed of . . . 
9000 words per second/per channel. Hence, the number of sampling points 
within the blade passage would be 400 for ft = 450 RPM. Because of, 
data blocksize limitations, this number n Q is reduced to 50 by artifi- 
cial 1 y multiplying the frequency of the signals 3ft by 8 with the ability 
of the analog tape recorder to record at 3-75 IPS and display at 30 IPS. 

A nine track digital tape is used to transfer the data to the 
computer which has a blocksize reading capacity of 12384 words at a time. 
This consideration limits the number of sets of 50 points equally spaced 
along one blade passage to 80. Thus the number of data available for 
determining the time average velocity or turbulence intensity at any of 
the 50 passage positions is eighty. 

Hence the blade to blade distribution of time averaged velocity is 
obtained through the averaging of the 80 consecutive sequences: 
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The fourth channel of the ADC unit is used as a control for the, 
read out of the data by the computer. 


3-3.3- Accuracy Criteria and Limitations 

A rough estimate of the accuracy of the blade to blade distribution 
of time averaged velocity is derived from reference 7* It is assumed 
that each instantaneous value within the blade passage is statistically 
independent. This assumption is supported by the fact that each value 
belongs to a physically different blade passage and that the integral 
scale T of the process could not be more than one period. The accuracy 
obtained is then given by (see Appendix 1): 


2 _ HI 2 

N 


(23) 


e is the difference between the actual average and the expected value, 
divided by the expected value. T is the turbulence level, that is the 
Root Mean Square (RMS) value of the fluctuating velocity divided by the 
expected value E{V} of the velocity. N is the number of samples 
averaged upon. 

The turbulence measurements described in the following section have 
pointed out a turbulence intensity of 10 to 15& for the relative 
velocity W, 75 to 80% in the radial direction (based on mean radial 
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velocity), and 15 to 20% in the tangential direction. The corresponding 
estimate for the precision e is, according to equation 23 and with 
N = 80, 

’ 1.12 to 1 . 7% for W . 

8.35 to 3% for V 

k. 5% for V z 

l. 67 to 2M for V Q 

Theoretically, these values may seem small and quite acceptable. How- 
ever, during the course of the experiment some qualitative limitations 
were noticed. These limitations are: 

1) It was pointed out in section 3-2.2 that the relationship (22) 
has to hold in order for the samples taken to be at a fixed location 
between the blades. Digitizing rate and analog tape recorder speeds are 
known to be accurate. However the inducer rpm can only be held fixed at 
± 1/10 of a rpm of the nominal rotative speed (^50 rpm). This could 
have led to a corresponding maximum shift of 2.3° in the 9 location of 
the samples over the entire averaging period. This represents 2 per 
cent of the blade spacing. Since the velocity gradients are not very 
large, a change in 1/10 of a RPM during the sampling time is unlikely to 
introduce any appreciable error. 

2) One of the very important limitations of this method lies in 
the fact that a three wire sensor can only detect the magnitudes of the 
flow along each coordinate direction without providing any indication as 
to the directions of the velocity. This is particularly true for the 
radial component of the velocity which is expected to reverse inside the 
passage near the blade tip. 


3) "Finally, the limitations of the ADC unit contribute to the 
accuracy of the method through its minimum increment while digitizing a 
fluctuating voltage. The maximum positive voltage accepted is 4. volts, 
divided in 128 windows, thus fractioning it in 0.031 volts increments 

which correspond to jump in velocity of 1.5 ft/sec for a fluid velocity 

* 

fluctuating around 20 ft/sec or 2.3 ft/sec for a fluid velocity fluctu- 
ating around 40 ft/sec. As can be seen from Figures (13), (14), (15), 
(16), the velocity profiles are not smooth even through the combined 
effect of turbulence and averaging. 

Turbulence Intensity Measurements 

The method of signal processing explained earlier is extended to 
turbulence intensity calculation. Since the instantaneous values of the 
three components of the velocity together with the averaged values are 
known for each sampling point within the blade passage, the root mean 
square of the fluctuating components of the velocity divided by the 
averaged value at each point are computed by the program. 

As in most turbulent cases, the probability density of the velocity 
is approximately Gaussian so that an expression for the accuracy of 
those results can be derived (see Appendix 1); 

E -\ff ^ 

E is an estimate of the accuracy of the measurements , giving 16% accuracy 
for N = 80 . 
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3-3. Blade to Blade Distribution of Time Averaged Velocities at 
The Exit . -.iMc • , 

The blade to blade variations of the three components of ’the abso- 
lute velocity at station 2a(Fig. ( 1 0) ) are obtained at seven radii; 

The values at four typical radii near hub (R = .582), mid radius (R = c 
.752), tip (R = .920) and halfway between mid radius and tip (R = .860) 
are shown plotted in Figures (13, 1^, 15, 16) respectively. In these 
figures the blade locations are identified on the basis of the measure- 
ment of turbulence intensities which are likely to be maximum at the 
center of the wake. 

1) Radial Velocity 

The blade to blade distribution of radial velocities derived from 
hot wire measurements are plotted in Figures (13a, lAa, 15a, 16a') for . 

R = .582, .752, .860, .960 respectively. The radial velocity distri- 
bution near the hub (R = . 582 ), plotted in Fig. (13a), r is almost 
linearly decreasing from pressure to suction surface except in a small 
region near the blade surface. The radial velocity outside the blade 
surface region can be attributed to blade blockage. At R = 0.752, the 
radial velocity distribution normally encountered inside the boundary 
layer of a rotating blade can be observed (Fig. (14a)). The radial velo- 
city is maximum very near the blade surfaces and is of appreciable mag- 
nitude at all other regions. The trend in radial velocity observed gives 
no indication as to the presence of a radial inward flow. 

At R = 0.860, the maximum radial velocity occurs near the suction 
surface. It should be emphasized here that the hot wire sensors can 
detect only the magnitude and not the direction of the radial velocity. 
Hence it is quite probable that a small region within the passage beyond 
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the pointrA: (Fig^.u^ 15a) ) has radial, inflow as compared to radial outflow 
at all other regions. The tuft probe mounted at the exit does reveal 
the. presence >of»:radi al i nward flow at this radius. Thus the peak 
observed at point B in Fig. (15a) might represent the magnitude of the 
maximum radial inward velocity. 

.'The radial velocity distribution plotted in Fig. (16a) for R = 0.920 
shows the presence of large radial inward, flow.. The maximum radial 
velocity region corresponds to the mid passage. In an earlier report 
(reference 5), a possible explanation for the existence of minimum rela- 
tive velocity inside the blade passage is given. This region corresponds 
to the location where the radial inward flow is maximum. 

•The blade boundary layer region in which the radial outflow occurs 
is'smal.l. This is confirmed by. rotating probe measurements to be 
described later, (see Fig. (17b)). 

Also apparent from Fig. (13a, 14a, 15a, 16a) is the fact that the 
magnitude of the radi.a.l velocity increases from hub to tip. It can be 
seen from Fig. (l^a) that the suction surface boundary layer is larger 
than that of the pressure surface. Also the radial outward., flows inside 
the pressure and. suction surface boundary layers decrease towards the 
tip (compare Figures (15a) and ( 1 6a ) . 

The radial velocities very near the suction and pressure surfaces 

derived from the measurement of the limiting streamline angle e (see 

w 

section 3 - 5 • 1 ) are shown plotted in Fig. (14a, 15a, 16a) and compared 
with hot wire measurements. The discrepancy observed is within experi- 
mental error. 


It is to be noticed here that the gradient in radial velocity 


av 

r 

r30 


is not as large as measured by Jabbari (3) on a rotating helical blade. 
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Since the turbulence level is high, the boundary layer wakes are likely 
to diffuse very fast as the flow proceeds downstream, thus considerably 
decreasing the velocity gradients in the tangential direction. 

Furthermore, the radial velocities plotted in Fig. (13a), (14a), 
(15a), ( 1 6a ) are of the same order of magnitude as the axial velocities. 

2) Absolute Tangential and Axial Velocity Distribution 

The absolute tangential velocity near the hub (R = .582), mid radius 
(R = .752), halfway between mid radius and tip (R = .860), and at the 
tip (R = .920) are shown plotted in Figures (13b), (14b) , (15b), ( 1 6b) 
respectively. The distribution observed is generally uniform with maxi- 
mum tangential velocity at the blade locations from hub to mid radius. 

At R = 0.86 and 0.92 (Fig. 15b, 16b) where there is a likelihood of 
radial inward flow outside the blade boundary layers, the tangential 
velocity is maximum wi th in the blade passage. Also shown plotted in 
Figures (13b), (14b), (15b), ( 1 6b ) are the axial velocity distributions. 
The blade to blade profiles are also surprisingly very uniform, with 
minimum axial velocities at the blade locations, as apparent in Fig. 
(14b) . ■ 

3) Relative Velocity 

The relative velocity distributions from blade to blade at the four 
radial positions reported here are shown plotted in Fig. (13b), (14b), 
(15b), (16b). From hub to mid radius the maximum relative velocity is 
located approximately at mid blade passage, whereas from mid radius to 
tip the maximum occurs near the blade locations. Also plotted on these 
Figures is the relative velocity distribution obtained from rotating 
probe measurements (see section 3.4). Although there is a discrepancy 
in magnitude between the two measurements especially near the tip, the 
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trend 6bs : erve'd’ v 'i £ E s iml lar. The reason for this discrepancy between the 
t^o'meas^u femed^s* J is given in section 3 • 3 • 5 . 

. ro'' js'C. i b _ < 

< ; H.*3 .3 . 5 • • Passage Averaged Values of Velocities at the Exit 

The blade to blade velocities discussed in section 3*3-^ are 

averaged over the blade passage. These are'plotted in Fig. (12a), (12b) 

andcompared with stationary and rotating probe measurements reported 

later. The good' agreement between the axial velocity profile at station 

2a derived from instantaneous velocity measurements and stationary probe 

(section 3-2) is evident from Fig. (12a). Also evident from this figure 

Js t: th'e "fact' that ‘the radial velocity measured is of the same order of 

magnitude as the axial velocity and increases from hub to tip. It is to 

be noted here that the average refers to I v [ and not v . 

' r r 

The magnitude of averaged radial velocities ( | v r | ) plotted in Fig. 
(T2a) confirms^ the earlier conclusion that the viscid and real fluid 
effects account for much of the flow characteristics at exit. Shown 
plotted in Fig. ( 1 2b) is the radial distribution of absolute tangential 
velocity at station 2a compared with stationary probe measurements. 
Although the trend of the two measurements is the same, with a large 
positive gradient near the tip, there is appreciable discrepancy between 
the two measurements. - Part of this discrepancy is due to the fact that 
the stationary probe reads the resultant of the tangential and radial 
velocity (when the radial velocities are small compared to tangential 
velocities) without providing any means of deriving them separately. 

This accounts for about half of the discrepancy between the two measure- 
ments. The remaining discrepancy must bei attr i buted to the experimental 
error. One possible source of error is due to the sensitivity of the 
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hot wire sensing V to the deviation from the cosineHato (section. 3.:3.- 1) . 
zr 

The corrected data derived from reference (13) lacks in accuracy 'for they! 
particular extreme case of 1/d = 700. '<■> ...’. i 

Also shown plotted in Fig.' (12b) -is the relative velocity at ' * 'vs. ; 
station 2a derived from instantaneous velocity measurements, compared to ' 
rotating probe measurements at station 2. The design value is also 
shown together with the corresponding rotating probe measurements at 
station 2 for the four bladed inducer. The discrepancy between 
instantaneous and rotating probe measurements can be partially explained 
by the flow redistribution that takes p.lace between station 2 and 2a. 

Nevertheless, an appreciable decrease in re,lat i ve . ve loci ty can be. : . 
noted between 3 and 4 bladed inducers at exit. This is consistent with 
a corresponding increase in static head. A large loss region is still 
present near the tip although it is somewhat reduced from that of a four 
bladed inducer. , ... 


3.3-6* Blade to Blade Distribution of Turbulence Intensity . 

The ; instantaneous and time averaged values of the three velocity,, 
components reported earlier are used to derive the turbulence components 
in all the three directions at various locations of the blade passage. 
The definitions for turbulence intensities used in this thesis are as 
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The^accurac/ of these turbulent intensities which are based on only 
eighty samp les i ls?g I ven in section 3 - 3 • 3 •, . 

The turbulent intensity components are shown plotted from. blade to 
blade at.R =.< i-752y v .860, .920 in Figures (13c and d, 14c and d, 15c and 
d?} 16c- and d) respectively. • 

The Radial Component T^ . , 

The radial component'of the turbulence 1 nteVis i ty^ shown plotted in 
Fig. J (13d), (lA'd) (15d) , (16d) for R= .582, .752, .860, .920 respec- 
tively is nearly unifornrfrom blade to blade. However, near the hub 
(Fig. (13d)) 1 the maximum is well marked at the blade locations and the 
turbulent intensity increases from suction to pressure surface, whereas 
near the'tip the distribution is more uniform (Fig. (16d)). The passage 
averaged turbulent intensity in the radial direction increases from 
about 70% near the hub to 75% at mid radius and tip. A slight maximum 
of 78% is found between mid radius and tip, Fig. (15d). This is consis- 
tent with the observation that there is large flow mixing due to the 
interaction between the pressure and ‘suction surface boundary layers 
somewhere between mid radius and tip. 

The Axial Compone'nt T^ 

V 

The axial component of the turbulence intensity is shown plotted in 
Fig. (13c), (14c), (15c), (16c) at the same radial • stations as above. 
Near the hub peaks in turbulent intensity are apparent at the blade 
locations, whereas near the tip it appears to be randomly distributed. 
Turbulence intensity in the axial direction appears to be approximately 


53 



\ 

1 
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constant from hub to tip around 40%, with a slight rnaximum of 43% in the 
boundary layer interaction region located between mid radius and tip. 

* S \ % 

\ ‘ ' 

The Tangential Component T Q 

■ £ . \ . 

Shown plotted in Fig. (13c), (14c), (15c), (16c) is the tangential 

H 

component of the turbulent, intensity. Much more uniformly distributed 
than the components in the radial and axial directions, the tangential, 
component is approximately constant from hub to tip around 20%, with a 
slight maximum at the blade locations. 

Turbulent Intensity T 
w 

The turbulence intensity of the relative flow is shown plotted In 

* 

Fig. (13d), (Hid), (15d), (16d) at the same four radial stations from' 
hub to tip. The distribution seems to be very uniform from blade to 
blade across the wake. The intensity is approximately constant from hub 
to tip around 14%. 

One of the major conclusions from these results is that the flow is 
general 1 y / very turbulent. The radial and axial components of the 
turbulence are maximum somewhere between the mid radius and tip. It was 
anticipated that the mixing and large radial motions due to boundary 
layer interactions at these locations of the blade passage would show 
a substantial change in flow properties from the rest of the blade 
passage. This is confirmed by this experiment. 
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3.^. Flow Survey Inside the Blade Passage Using Rotating Probes 

f 

To find the flow properties inside tb.e blade passage the flow is 
surveyed at station 2 (Fig. (10)) using a rotating pitot tube and a 
pressure transfer device described in sections 3.1.2, 3-1.^. The 
stagnation head (\^j) R of the relative flow from blade to blade is deter- 
mined at ten radial locations at station 2. The static pressures ij; at 
corresponding radii are taken from holes drilled on the pressure and 

A 

suction surfaces. These static pressures are assumed to, vary linearly 

from pressure to suction surface. The relative veloci ties- are then 

derived from (>Jj t ) and . ‘ 

I K S 

, t 

The blade to blade distributions of the stagnation head of the 

\ 

relative flew and the static pressure coefficient at lOvradial stations 
from hub to tip, at station 2 (Fig. (10)), are shown in Fig. ( 1 7a- j ) . 
Also shown plotted on these figures is the corresponding Relative 
velocity. 

• The usual flow turning effect in addition to the boundary layer 
growth on the blade can be seen from the relative velocity plots' for 
R = 0.506 through R = 0.720 (Figures 17f - j). The boundary layer 
growth on the suction surface is generally higher than that of the' pres- 
sure surface. However it seems that the extent of the pressure and' 
suction surface boundary layers have diminished considerably compared to 
the four bladed inducer (Ref. 5), Figures (1 7a) , ( 1 7d) , ( 1 7j ) • From’hub 
to about mid radius Fig. (17f - j)» the blade to blade distribution' of - ' 
relative stagnation head is uniformly distributed and rather constant, 
whereas from mid radius to tip it decreases from pressure to suction 
surface Fig. (17a - e) . 
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The relative velocity increases from pressure surface to suction 
surface from hub to mid radius due-to the usual flow turning effect • „= 

(Fig. 1 7f - j), where as from mid radius to tip a minimum is noticeable 
within the passage located nea,rer to the suction surface (Fig. 17a - e) . 
Also, it can be seen from Fig. ( 1 7d) that the trend of the relative 
velocity is unusual" in that* it decreases from pressure to suction sur- 
face. This distribution of relative velocity is unconventional, and is r 
probably caused by the i'nteract-i on 'effect between the blade boundary ... . 
layers which generates’^- large radial inward flow within the blade pas- 
sage. The variations in radial velocities from blade to blade are thus 
large, causing the type of profile observed in Fig. (17a - d) . This 
effect is most noti ceable at radius R = .808,Fig. (17d), where the radial 
inward flow is ITkely to be maximum wi thi n the blade passage. Similar 
conclusions are' reached; from instantaneous velocity measurements and : the 
relative velocities obtained by both methods are shown compared. in 
Figures (17b v ,17e, and 17h). it is to be observed that radial velocities 
and 'turbulent-' intensi ties (T , T ) are maximum here (See Figures 15a, . 

* K Z 

15c). Although the trend observed is the same, the discrepancy between' 
the two measurements has been explained in section 3.3-5, and is partly 
due to experimental error. 

J The three bladed and four bladed inducer results are compared in 
Figures* ( 1 7a , d, j). It can be observed from these figures that while, 
the static pressure coefficients have increased from three to four : 
b-laded inducer, the stagnation head of the relative flow has nearly .s 
remained the same. Consequently, a decrease in relative veloci ty . is 
observed i n the case of three bladed i nducer (See Figures ?( 1 7a, : d; j)'. ' 
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Also, the head and relative velocity profiles are more uniform than in 
the case of the four bladed i.nducer from^hub ,to mid radius. 

3- 5- Flow Visualization Method 

As. explained in reference a knowledge , of the limiting streamline 
angles on the blade surface provides very useful information on the 
nature and magnitude of the radial flows and the direction of the wall 
shear stress. With a view to obtain this information for a, future 
theoret i cal. sol ut (on of the flow, an attempt is made to derive these 
angles by ammonia streak technique and smoke produced by titanium 
tetrachloride. 

1) Limiting Streamline Angle on the Blades. 

The method of deriving the limiting streamline angle on the blades 
at the. leading edge- has been described in reference 3- Absolution of 
titanium tetrachloride and carbon tetrachloride is prepared. When 
exposed to the airflow, this solution generates a white, dense smoke 
capable of being photographed. The mixture Is brushed on the, surface of 
the blade and photographs are taken while the blade is rotating. The 
deflection: of the smoke filament from the tangential direction gives an 
indication of e^. The results of this measurement at the leading edge 
are shown plotted in Fig. (18) and compared with four bladed inducer 
result's. e w is seen to be decreasing 1 inearly- f rom hub to tip, where its 
value reaches nearly zero. Surprisingly at the leading edge is lower 
than in the case of the four bladed inducer. 

The second method to derive the limiting streamline angle at exit 
is the ammonia str.eak filament method. The static pressure holes 
drilled on the surface of the blade (see section 3- 1.1) are used in 
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combination with the ammonia transfer device. A small amount of ammonia 


gas at low velocities is fed into the tubes at various tangential and 
radial locations of the blade surface. A sheet of ozalid paper sensi- 
tive to ammonia is pasted along the edges of the static holes and traces 
of ammonia are recorded on the ozalid paper while the inducer is in 
rotation. The results of the measurement of at trailing edge on . 
suction and pressure surface are shown in Fig'. (18) and compared with’ 

the four bladed inducer results, e on the suction surface is seen to 

w 

decrease from 50 ° at mid radius to zero at tip near the trailing edge. 

Its value is generally higher than the corresponding e for the four 

bladed inducer, where - as e is lower on the pressure surface.- Further- 

w 

more, the magnitudes of are generally higher than those observed on la 
single helical blade, (reference 3). These results suggest, that-' radial - 
velocities encountered in inducer passages, where cons i c de-rab le velocity 
and pressure gradients exist in all directions, are generally higher- ■ 
than those observed on a single helical blade. 

2) Tuft Grid at Exit % . 

The maximum deviation of the streamlines from the- cy 1 h ndr-J caile.no « 5 * 
surface (e ) , outside the blade boundary layer regions, and obtained 
from a tuft grid mounted at exit of the inducer is shown plotted in ' * 
Fig. (18). Although there are appreciable radial inward flows from mid 
radius to tip, their magnitudes are generally lower than those of the 
four bladed inducer. This confirms the earlier observation that the- 2 
blade boundary layer's interaction effects can be reduced by decreas.ing; 
the solidity of the inducer. 25 ' >n-' : - o J ■ 
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4. DISCUSSION 


Three Bladed Inducer Performance 

Axial Velocity 

The outlet radial distribution of axial velocity derived from 
stationary probe measurements shown plotted in Fig. (19b) at several 
axial stations indicate that the backflow region near the hub originates 
at the trailing edge and extends downstream. The flow separation zone 
reaches 50% of the annulus region far down stream of the inducer. A 
large flow redistribution is seen to occur between -stations 3 and 4 , 
probably due to large axial gradients in radial velocity. The axial 
velocity perturbation at exit obtained from the blade blockage analysis 
of section 2.1, plotted in Fig. 19b, seems to indicate that blade blockage 
effects are not negligible in the case of the Penn State inducer. 
Furthermore, the qualitative trend of the distribution is predicted at 
stations 3 and 4, thus accounting for the part of the change in axial 
velocity downstream. The agreement between instantaneous velocity and 
stationary probe measurements for the axial velocity Is good (Figure 12a). 

The comparison between measured axial velocity at station 4 and the 
predicted values derived from the exact inviscid analysis described in 
section 2.2.2 shows that the three dimensional inviscid effects on axial 
velocity are not small. However, the three dimensional inviscid solution 
does not predict the back flow near the hub and the very steep gradient 
near the tip. This shows that an accurate prediction of the flow at the 
exit of inducers has to include viscid effects. In an earlier report 
(references 4, 5) it is argued that all the deviations in axial velocity 
are due to viscid effects. The present analyses tend to indicate that 


59 


this is not quite true and that three dimensional inviscid effects 
account for some of the deviations 1 ; The blade to blade variation of 
axial velocity is very uniform at all radii from hub to tip, with a 
slight minimum at the blade location, Fig. (13b). The axial component 
of the turbulence intensity from blade to blade at exit and at several ' 
radii is a 1 so surpr i s i ng 1y uni form throughout the blade passage with 
value around ^ 0 %. 

•Radial Velocity v 

The averaged radial 1 ‘ve loci ty ( | V r | ) at station 2a increases from 
hubto tip, Fig. - (12a) : : The radial velocities, derived from three 
dimensional hot wi re‘ measurements', are found to be"- of the same order of 
magnitude as the axial velocity in the entire passage. The departure 
between the predi cted- (Fig . (**)) and the measured radial velocities at 
outlet is very large. A major reason for' the existence n of the radial 
velocity in the passage is due to viscid and real fluid effects. The 
blade to blade distribution of radial velocity between mid radius and 
tip, Fig. (15a), (16a), is found to be characteristic of three 
dimensional boundary layer interaction effects. The flow vi sual i zati on" 
indicates the presence of large radial outward flow within the blade 
boundary layers and radial inward flow within the blade passage from 
mid radius to tip. These measurements confirm the reliability of the 
' radial" ve loci ty profiles derived from hot wire measurements (see Fig. 
(15a), ( 1 6a ) ) . These results show that the radial velocities encountered 
in inducers are of the same order of magnitude as the axial velocity, ' 
at least near the exit of the blade passage. ’’ '• 
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Absolute langentiai velocity 

The radial distribution of absolute tangential velocity indicates a . 
steep rise in the values near the tip (Fig. (19c)). Furthermore, there 
is considerable change in the tangential velocity as the flow proceeds 
downstream. Fig. {k) shows the comparison between the predicted tan-. • 
gential velocity at exit using the three dimensional analysis of section 
2.2.2, and the experimental value at station 3. The predictions are 
good from hub to mid radius while the discrepancy is large near the tip. 
Hence it is obvious that accurate prediction of the tangential velocity 
predicted by the approximate viscid analysis of section 2.3 agrees well 
with the experiment (Fig. (8)), thus showing the validity of the pres- 
sure loss coefficient derived in reference (4) and plotted in Fig. (7). 
There is some discrepancy between the tangential velocities derived from 
hot wire measurements and the stationary probe, Fig. (12b). Part of it 
is due to the fact that stationary probes sense both the radial and 
tangential velocity and the remaining could be attributed to the experi- 
mental error, since the correction factor for the deviation of the cosine 
law. used. in hot wire measurements is not known very accurately for this 
particular case. The blade to blade distribution of absolute tangential 
velocity is found to be uniform throughout the blade passage at exit, 

Fig. (13b), ( 1 Ab) , (15b), (16b), with maximum values at the blade ■ 
location from hub to mid radius. The maximum tangential velocity occurs 
outside the blade boundary layer near the tip region. The tangential 
component of the turbulent intensity has a maximum value of 22% in the 
boundary layer interaction zone which is located between mid radius and 
tip- 
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Relative Veloci ty 

The relative velocity increases from hub to about mid radius and is 
very close to the design values (Fig. (12b)). From mid radius to tip, 
the departure between the design and measured values is very large thus 
indicating a major flow departure and large mixing losses in this region. 
The discrepancy between the rotating probe and hot wire measurements for 
the relative velocity is due to the differences in the methods explained 
earlier. It should also be noted that the rotating probe measurements 
are taken within the passage at station 2, while the hot wire measure- 
ments at station 2a are taken outside the blade passage. Thus, part of 
the discrepancy is more apparent than real due to the flow redistribution 
The blade to blade variation of relative velocity, shown plotted in Fig. 
(13b), ( 1 -4b ) , (15b), ( 1 6b ) for four radial stations at exit, is nearly 
uniform indicating the fast decay of the blade wakes. From hub to mid 
radius, the maximum value for the relative velocity seems to occur near 
the blade locations. 

The turbulence intensity of the relative flow shows a uniform trend 
near the tip, whereas maximum values at the blade location are notice- 
able near the hub, (Fig. 13d)). The average value of the turbulence 

intensity T is found to be about 13%. 
w 

The *.'re 1 iab i 1 i ty of the blade to blade distribution of the relative 
velocity is confirmed by the comparison between hot wire and rotating 
probe measurements at three radial stations, Fig. (17b), (17e), ( 7 7h) . 
Although the discrepancy in magnitude noted earlier is there, the trend 
observed is the same. These measurements have provided some valuable 
information for the theoretical prediction of the flow field including 
viscous effects. 
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Comparison of Three and Four Bladed Inducer Performance 

From the results reported in this thesis and reference 5, it is 
found that the three bladed inducer has appreciably better performance 
than the similar four bladed inducer tested at the same flow coefficient. 
In particular, the following conclusions can be drawn from the investi- 
gation of the flow in three bladed inducer. 

1) The stagnation and static head rises are appreciably higher 
than those of a four bladed inducer at all radii. Even though the steep 
rise for the head coefficients towards the tip is still present, their 
gradients are somewhat reduced. Further improvements may be possibly by 
boundary layer control on the blades. The increase in the values of ^ 
in the case of a three bladed inducer from those of the four bladed 

,J -f 

inducer is of the order of 15 to 30 per cent. 

2) The axial velocity profiles show a marked deterioration from 
that of a four bladed inducer. There is an increase in the extent of 
backflow near the hub compared to that of a four bladed inducer. The 
axial velocities near the blade tip are found to be generally higher 
than those of a four bladed inducer. 

3) There is no appreciable difference in absolute tangential 
velocities measured at the exit of the two inducers. This indicates an 
improvement in hydrodynamic efficiency which is defined as the ratio of 
the actual head to the Euler or ideal head rise. 

**) There is considerable improvement in hub wall and annulus wall 
static pressure di s t r ibut ion. The negative pressures observed near the 
leading edge of a four bladed inducer are appreciably reduced in the case 
of a three bladed inducer, Fig. (20). Furthermore, the negative radial 


63 



pressure gradient measured In the case. of a four bladed inducer at some 
axial stations has completely vanished in the present case. 

•5) The; relati ve, veloci ties., at the.., exit are gene rally-, .lower than 
the four bladed inducer, especially near -the tip. -Th;is is consistent 
wl th the -i ncrease in. s tati c pressure observed in. the- -case of three' 
bladed , inducer-. ...The flow. mixing losses: near the tip appear to have 
decreased.*, - . . :• 

. 6)^ . There is substantial ..improvement i n the uniformity of the blade’ 
to blade distribution of relative velocities inside the^blade passage.- ' 
The. extent of- the blade boundary layers seem, to' have appreciably 
dimi ni shed. .... - ; -•• • ■; 

-.7) .The ampl i tude of the radial inward flow wi th i n . the b.lade 
boundary layer interaction region seemsto have been reduced considerably, 
in the case of the, three bladed Inducer. . . 

Further . improvement, of the flow, characteristics, .could- probably -be 
achieved by boundary layer-control, on th.e^blade and by a reduction in .. 
blades thickness. An improved analysis for the prediction of the u >flow 
field- at exit should be based on a better knowledge, o ; f- the radial 
velocities .inside the-blade passage. . -Re-1 iable information on these . 
radi al veloci t ies could come- from the study of the three-dimensional 
boundary- layers in a rotating helical channel, currently under Investi- 
gation at. the Department of Aerospace Engineering of The Pennsylvania - * 
State University. ■ * 
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5. CONCLUSIONS 


The following conclusions may be drawn from the theoretical and 
experimental' investigations reported in this thesis. 

1) The perturbation of the meridional flow caused by the radially ' 
varying blade blockage is not negligible for the Penn State inducers.- 
The analysis carried out in section 2.1 predicts perturbation of up to 
10* in axial velocity at the exit. The axial velocity profile derived 
from this analysis agrees qualitatively with the experiment. 

2) The program developed by P. Cooper and H. Bosch for the exact 
solution of the inviscid flow field has been somewhat improved to satisfy 
the Kutta Joukowski condition at the trailing edge and to include the 
blade blockage effect. The axial velocity profile at the exit derived 
from this analysis agrees qualitatively with the experiment. The 
tangential velocities aiso agree qualitatively with experiment, except 
near the tip where viscous and turbulent effects are dominant. 

Thismethod of solution could be extended to include real fluid 
effects by including the dominant shear stress terms in the equations 
of motion. An ini tial attempt was made in this direction by the use of 
the boundary- layer data and the skin friction coefficient derived by 
Jabbari in reference 3. In view of the large number of grid points 
required and the considerable increase in computer time, no attempt was 
made to run this program. 

3) The approximate viscid analysis of the flow carried out in this 
thesis is based on the assumption of fully developed turbulent flow in- 
side the passages and takes into account the frictional losses. The 
friction loss coefficient applicable to inducers operating in the range 
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of flow coefficients $ = 0.065 to 0.2 derived empirically -by.'is" 

Lakshmi narayana in reference 5 is used. This friction factor is’ : found 
to increase exponentially towards the tip. The values ofctthe. relative . 
and absolute tangential velocities predicted from thi's analysis agree 
closely with the experimental values. However , -further refinements in 
the radial velocity estimate are necessary for the accurate prediction 
of the axial velocities at the exit of the inducer. 

4) A method measuring the three velocity components and the corres- 
ponding turbulent intensities is developed. Th i s' procedure wh i ch 
utilizes three hot wires located in the coordinate directions has led to 
valuable information on the blade to blade variation of all the velocity 
and turbulence intensity components. More precise information on the 
deviation from the cosine law for the hot wire probes should be. obtained 
since this is one of the many sources of error involved in this method. 

5) The radial velocities obtained from the hot wire measurements 
are found to be of the same order of magnitude as the axial velocities 
throughout the flow passage. Near the hub, the radial velocities are 
found to be nearly uniform across- the passage and this is presumably 
caused by the blade blockage. At other radial 1 oc’at ions',' the radial 
velocity distribution across the passage is of the type norma My 
encountered in rotating boundary layers. 

6) Almost identical axial velocity profiles at the exit are 
derived by both hot wire and stationary probe measurements. There is 
considerable decay in axial velocity profiles in the case of the three, 
bladed inducer. ’ 

The blade to blade variation of axial veloci ties-measured at the 
exit by means of hot wire probes is found to be nearly uniform. 
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7) Appre.ci able discrepancy is observed in the absolute tangential 
velocities at 'the exi t measured by stationary and hot wire probes. Steep 
rises in tangential velocity are present. ,i.n the case of the three bladed 
inducer even though the gradients are somewhat lower than those observed 
for the four bladed inducer. 

The blade to blade variation of tangential velocities measured by 
hot wire probes indicate peaks near the blade locations from hub to mid 
radius where as at other radii the peaks occur inside the passage. The 
blade to blade variation of tangential velocity is found to be generally 
uniform across the passage. 

8) There is slight discrepancy in the blade to blade variation of 
relative velocities measured by the hot wire probe and the rotating 
pitot tube. However, the trend in their distribution is the same. 

The passage averaged relative velocities measured by both methods 
show a considerable departure from the design values from mid radius to 
tip. These averaged relative velocities are somewhat lower than those 
of the four bladed inducer at these radii. 

The blade to blade distribution of the relative velocities measured 
by the hot wires is found to be nearly uniform. This is probably due to 
considerable wake diffusion that takes place between the trailing edge 
and the hot wire measuring station. 

9) The stagnation and static head rise of the absolute flow are 
found to be appreciably higher than those of the four bladed inducer at 
all radii. The steep rise in stagnation head coefficient is still 
present in the case of three bladed inducer, even though the gradients 
observed are somewhat reduced. 
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10) There is appreciable improvement in the hub and wall static 
pressure distributions. The negative pressures observed near the leading 

.. • -l i . _■ : -:t-l 

edge of a four bladed inducer are appreciably reduced in the case of a 
three bladed inducer. 

'll) The turbulence intensity components in the radial, axial, 
tangential and relative directions (based on local mean velocities) are 
found to be about 75%, 45%, 20% and 15% respectively. Although a maxi- 
mum in radial and axial turbulence intensity is observed in the boundary 
layer interaction region at R = 0.860, the turbulence intensities are 

• i . ■ , ; 

found to be rather uniform in the entire flow passage. These measure- 
ments reveal the highly turbulent nature of the flow in inducers and 
should be extended to the determination of Reynolds stresses. 

12) The measurement of the relative flow inside the blade passage 
obtained by means of a rotating pitot tube also shows the two different 
flow patterns across the blade passage; from hub to about mid radius the 
blade to blade distribution of relative velocity and stagnation head of 
the relative flow show the usual trend characteristic of flow turning. 
From mid radius to tip the interaction effect between blade boundary 
layers becomes apparent and a wake type of profile is seen. This measure- 
ment has revealed an appreciable decrease in blade boundary layers as 
compared to the four blade inducer. 

13) The radial inward flows observed by means of tufts mounted at 
the inducer exit indicate that their magnitudes are generally lower than 
those observed at the exit of the four bladed inducer. 

I**) The limiting streamline angles measured on the blade surfaces 
of the three bladed inducer follow the same trend as four bladed inducer 
even though their magnitudes are generally higher. 
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15) Some improvement in exit flow characteristics of axial flow 

(low 'it.'- ■ - : 

inducers has been achieved by reducing the solidity. The three bladed 
i nducer studied in this thesis presents higher and more uniform head 
characteristics at exit than the four bladed inducer. The influence of 
viscosity and the interaction effect between the blade boundary . layers 
have considerably decreased. However a deterioration in exit axial 
velocity is observed, together with a larger flow separation zone 
extending downstream of the trailing edge. Further improvement of . the 

* ’ C * * - * • » ' 

inducer should include blade boundary layer control and possibly a 
blade thickness reduction. 



APPENDIX I 


Derivation of Accuracy Criteria for Instantaneous Veloci ty~ Measurements 

(Equations 23 and 24) 


1) Accuracy of the blade to blade averaged velocities (equation 23) 
If N successive samples are taken at any particular location of 
the blade passage (9), the value of the velocities at these samples 

are: V V v 3 V 

Taking the average: 


“ 1 ? 

V = rr E, V 

n N n=l n 


with the assumption that those n values are statistically independent, 
i t can be wri tten: 


N 


V n “ V 


N 


rr E. (V ) - V = rr E. (v - v) 
N n=l n N n=l n 


v being the expected value of the sequence v^ , . . 

(or the value to which the sequence converges for n “ °°) 
The standard deviation is given by: 


N 

E 


E {v - v) 2 } = — , E E (v - v) (v. - v)} 

n N 2 n , 1 = 1 n 1 


the independence of the measurements gives, 

E {(v - v)> E {(v. - v)} =0 

n 1 


2 i N 

and E { (v - v) } = — =■ E. E { ( v - v) } 

n n= I n 
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Dividing by v 


.2 

£ , , . o 


N 

E 

ri=l 


- v 


E ) 


hence 

e 2 * (a/x) 2 or e 2 = ^ (T) 2 

Note: the notation E{x} accounts for: expectation of the value of x. 


2) Accuracy of the turbulent intensities measurements 
Since the average is taken over N samples, the square of .the 
relative error for the fluctuations is 


E{ (v - \T) 2 } 
n n 

2 

v 


or 


N 

E 

n= 1 


that is, if 


E{(v p - v) 2 } 


E{ (v -v~) 2 
n n 

2 


} 

= 2 for a Gaussian process, 


E{(v-vJ 2 } , 

0 n 2 

2 = N 

v 

And the accuracy with which the turbulent intensities are obtained is: 
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APPENDIX 2 


Fortran Listing of the Program Used for Deriving instantaneous Velocities 

From Hot Wire Measurements 

// tape=alain, m it 

// FND flF TAPE CARDS 

// EXEC FGCLC,,PARM.ORJECT = MAP 

//source. input on * 

I 

DFRIVATI ON OF THE THREE COMPONENTS OE THE VELOCITY IN THE WAKE 
OF THF PENN STATE THREE BLADE INDUCER AT STATION 4 

■ f . r 

PE Al.*8 Q( ?,.101 ) , 

R F A i. 7n(4,>000) , „ . 

. DIMENSION Ai(51,80),Clf51,80),BI(51,80>,01(51,80),UT(102),UZ<102)» 
KIR (10?)', W( 102 ) , TUT ( 10?) ,TUR(.l02) ,TUW( 102 ) * TUZ ( 10?) , WT ( 5 1 ) , CVf 102 ) 

1 ,CZ ( 1 0? ) ,CR ( 10? ) ,CW( 102) , 

IWTFGFR*? I Al.f 2000) , I RL( 2000) , !CL( 2000) . " * 

E0U1VAI.FMCE ( L ( 1 ) , I NT ) , ( CC ( 1 ) ,01 1, 1 ) ) . ’ 

LOO I C Ai.* 1 CC('l f 63R4) ,1.(4) " 

READ 1002,K,Y1 ,Y2,Y3,Hl,H2t H3,T1,T2,T3 
PRINT )007,K',Y1,Y2,Y3,H1,H2,H3,T1,T2,T3 
M=1 • ‘ ' 1 

• N=1 ’ ' 

C **THE DATA CORRESPONDING TO EACH RADIUS .IS READ** 

98 CALL GET (91, CC, 16384, KEY) 

IF (KFY) 300,101,301 

301 PRINT 1004 • ' 

GO TO 98 
300 PRINT 1003 
N='N+1 

I F (N-4) 98,98,99 
101 IF (N-ll ) 100, 100,99 
100 I NT = 0 ... , 

R F AD 1 008 , R 

IR=l . . 

I C = l . ,. 

C **THE DATA is PUT IN TWO DIMENSIONAL ARRAYS** 

DO 200 1=1,16320,4 

L ( 4 ) =CC ( I ) ‘ 

A.l (IC»IR) = 1 28—1 NT 
L ( 4 ) =CC < I + 1 ) 

R1 ( IC, IR)=INT-128 
L( 4) =CC( 1+2) 

Cl ( IC, IR)=INT-128 

L( 4) =CC( 1+3) - 

ni ( IC, IR) =INT-128 ’ ! 

1 F ( IC.LT.51) GO TO 200 • = ’ 

IC = 0 

I R= TR + 1 I , 

200 IC=IC+1 1 

PRINT 1021 
M = M+l 

PRINT 100? , M , ■, . .... ■ • 

RUP=0. 0 , | ' 

RMM = 0.0 '■ , 

Rl 10 = 0. 0 

C ** BEGINNING OF THE AVERAGING PROCEDURE** ' .. 

DO 3 1=1,51 

PRINT 1001,(I,A1(I,1),B1(I,1),C1(I,1),D1(I,1)) 

SUP=0.0 
SUM=0.0 
SU0=0,0 
SI IR = 0 *0 
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4 


3 


C 


5 


30 


C 


S I IM = SI IM+ A 1 ( 1 
SUP = SUP + R 1 ( I 
siiR = suR+m (,! 

R||P = RIIP + A1 ( ! 
RIIM = RIIM+ftl ( I 
Rlin = RIIO + C 1 ( I 

s ho= si io-i-c i < i 

AV1 = SIIM/R0. 0 
flV?=SUP/fiO.O 
AV3=SU0/R0. 0 


.1 ) 

J) 

J)l 

J) 

J ) 

J) 

J) 


i 


AV4 = /S4.0 

TAVl = ( AVl/AV4)*Yl' ! 

TAV?=( AV2/AV4)*Y2 
TAV3=('AV3/AV4)*Y3 

R 1 = I ( (TAV1*TAV1-Hl )*( TAV1*TAV1-H1 ) )/( Tl*Tl) ) 
R2=( ( (TAV2*TAV2-H2)*ITAV2*TAV2-H2) )/(T2*T2) ) 
R 3 = ( ( (TAV3*TAV3-H3)*( T AV3*T AV3-H3 ) )/<T3*T3> ) 
Xl=Rl*Rl/< 2.0*0. 9fl) , 

X2=R2*R2/( 2.0*0.98) 


X3 = R3*R3/< 2.0*0.<5fl) 

I IT ( I ) = ( SORT! ARSI X1-X2 + X3 ) > ) / 7 1 . IS 
UZ. (I) = < SORT! A RS ( 0. 9#»*Xl + X2-X3) ) 1/71.6 
HR { l ) = ( SORT ( ARSI X2 + X3^0.9ft-X 1 l ) ) ) / 7 1 . <S 
WT< I >=R-UT( I ) 

W( I ) = SORT ( UZ ( 1 )*UZ< 1 ) +UR ( I )*UR( 1J+WTI I)*WT< !) ) 


RUP = RIIP./40R0.0 

RIIM = RIIM/40RO.O 

RIIO = RUO/40«n.0 

PRINT 1001 , I ,RUP, RUM, RUO 

PRINT 1006 

PRINT 1 009 , R 

PRINT 1000 

**TIMF AVFRAGFO VELOCITIES** 
nn s i=i, si 

PRINT 1001 , ( 1 ,UZ( I ) ,UR(I) ,IIT(I) ,WT( I ) ,W(I ) ) 

SUM=0. 0 

SUP=0.0 

suo=o.o 
suR=o.n 
sus=o.o 
nn Rn i=i, si 
SUM=SUM + I)T ( I ) 

SUP = SUP-*-* I Z I I ) . ' • 1 

S • 10= SI 10+IIR ( I ) 

SUR=SUR+WT ( I ) 

sus=sus+wi I } 


SIIM = SUM/S l .0 
SUP=SUP/51 . 0 
SUO=SUO/Sl .0 
SIIR = SUR/S1.0 
SMS=SUS/Sl .0 
PRINT 1020 
PRINT 1000 

PRINT 1001 , K, SUP, SUO, SUM, SUR,SUS 

**TURRIILFNCE INTENSITY CALCULATION** 

00 21 1=1,51 

SUP = 0.0 

SUM = 0.0 

SUQ ■ 0.0 
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SDR = 0.0 

nn 22 . 1 = 1, 80 

AVI = A 1 { I , J) 

AV2 = Bl( I, J) 

AV3 = Cl( I, J) 

AV4=64.0 

TAV1 = ( AV1/AV4KY1 
TAV2=( AV2/AV4)*Y2 
TAV3=( AV3/AV4)*Y3 

Rl=( ( (TAV1*TAV1-H1 )*( TAV1*TAV1-H1 ) ) /( Tl*Tl) ) 
R2=( ( ( TAV2*TAV2-H2)*( TAV2*TAV2-H2) 1 / < T2*T2 1 ) 
R3=< ( (TAV3*TAV3-H3)*(TAV3#TAV3-H3) )/<T3*T3)l 
Xl=Rl*Rl/(2.0*0.9ft) 

X2=R2*R2/(2.0*0.9ft) 

X3=R3*R3/( 2.0*0.96) 

A = (SORTIABSIX1-X2+X3) 1 1/71.6 
B = <S0RT(ABS(X1*0.96+X2-X3) 1 1/71.6 
C = ( S0RT( ABS< X3*0.96+X2-X1 ) ) 1 /71 .6 
n=R - a 

F = SORT ( RSB+C^C+D^D ) 

Slip = SDP + ( A-UT ( I ) ) * ( A-UT ( I ) ) 

SUO = SDO + (B-UZU ) )*( B-UZ( I ) ) 

SUR= SUR+ ( C — IJR ( I ) ) * ( C-UR ( I ) ) 

22 SOM = SUM + ( E-W (!))*( E-W( I ) ) 

TDR ( I ) = { ( SORT) SDR /BO. 0) )/DR< I ) 1*100.0 
TDZ ( I ) = ( ( SORT ( SD0/80.0) ) /DZ ( I ) ) *100.0 
TDK I 1 = ( ( SORT < SDP/80.0) )/DT< I ) ) *100.0 
21 TDW( I )=( ( SORT! SDM/ 80.0) 1 /W| I) 1*100.0 
PRINT 1012 
PRINT 1013 
00 1ft 1=1,51 

lft PRINT 1001 , ( I » TDT ( I ) » TU-Z ( 1 ) » TDR ( I ) , TDW( I ) ) 

00 17 1=52,101 

J= I -5 1 

DZ( I )=DZ( J) 

DR ( I ) =DR ( J ) 

DT ( I ) =DT ( J ) 

W(I)=W(J) 

TDZ ( I )=TDZ( J) 

TDT ( I ) =TDT ( J ) 

TDR ( I ) =TDR( J) 

TDW ( I ) =TOW ( J ) 

17 CONTINDE 

DTI 101 )=0.0 
TDZ ( 1011=0.0 
W( 101)=0.02 
TDR ( 1011=0.0 

C ssGRAPHIC ROOTINES** 

IF (M-2) 98,305,309 

309 IF (,(M-5)*(M-6)*(M-7)*(M-B) 1 98,306,98 
305 CALL INI TO I I AL, I BL, ICL, ZD, 2000) 

30ft 00 ft 1=1,101 

0(1,11=1 

ft 0(2,1) =l)R ( I 1 

CALL GRAPHOC $',' S', 0,101,0) 
no 11 1 = 1,101 

11 0(2, I 1 =W ( I ) 

CALL GRAPHOI' $',0,101,01 

no 9 1=1,101 
9 0(2,1 1 =l)T ( I 1 
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**. 0 , 101 , 0 ) 


**, 0 , 101 , 0 ) 


CALL SAMFO 

CALL GRAPHQ ( * *',* **,0,101,0) 

00 R 1=1,101 
R 0(7.1 )=UZ( ! ) 
cai . i . samfo 

CALL GRAPHO | • $*,' **,0,101,0) 

00 13 1=1,101 

13 0(7.1 )=TIIZ( I ) 

CALL GRAPHO) ' **,* *', 0,101,0) 

00 17 1 = 1.101 

1? 0(2.1 ) = TOT ( I ) 

CALL SAMFO 
CALL GRAPHO(* 

00 14 1=1.101 

14 0(7. I ) = TOR ( I ) 

CALL GRAPHO! * 

no i5 1=1,101 

1 5 0( 7, I )=TIIW( I ) 

CALI. SAMFO 

CALL GRAPHOl • *',' **,0,101,0) 

1000 FORMAT < //2X, 1HI , 10X.2HI/Z, 10 X , 2H0R , 1 OX , 2HUT , 10X.2HWT, 10X, 1HW/) 

1001 FORMAT (' ' , I 5, 1P4E13.6) 

100? FORMAT ( 1 2.3F6.3.6E10. 3) 

1003 FORMAT (' '.24HTAPF FRR OR -PROGR AM STOPPEO) 

1004 FORMAT (• '.IIHFMO OF FILE) 

1004 FORMAT (ROH EXIT VELOCITIES NON D I MENS I ON AL I Z ED WITH RESPECT TO T1 
IPSPFEn AT RAOIDS LOCATION) 

1007 FORMAT (* ■ , I2,3F4.2,6E10.3) 

l OOR FORMAT (FI0.4) 

100Q FORMAT ( ' ' ,F 10.4) 

1010 FORMAT!* • , 15, 1P4E13.6) 

1011 FORMAT ( //2X, 1HI ,IOX,2HA1,10X,2HR1, 10X.2HC1, 10X.2HD1/) 

1012 FORMAT!' ' , 3RHTUR SOLENT INTENSITIES AT EACH LOCATION) 

1013 FORMAT! //2X, 1HI , 10X.2HIIT, 10X.2HIIZ, 10X.2H0R, 10X, 1HW/) 

1070 FORMAT (' '.26HRMS VALDES AT THAT STATION) 

1071 FORMAT (• 1 


99 


]*«•«**««**«*« ) 
GO TO OR 
STOP 
FNO 


/* 

//0ATA.FTR1F001 no UNI T« ( 2400, , DEFER ), VOL-SER- ALA IN, X 

// LARFL»< l.NU ,nCR»(RECFM=F,BLKSI Z E » 1 G384 , RUFNO- 1 ) 

//nATA. input no * 

7 0,999 1.350 0.999 .R20E+00 .730E+00 .760E+00 .177E 00 .165E 00 ,179fc 00 

. 5R 7 
,43R 
.49? 

.754 
• R09 
. R44 
.970 
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Figure 1: Photographs of Four and Three Bladed Inducers 


78 












view into inlet 











second run,exit direction free 
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suction surface 
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FigureSb; Blade Pressure Distribution at Mid-Radius 


pressure surface 
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F igure 5c: Blade Pressure Distr i but ion at TIP 




85 









from Approximate Viscid Analysis 
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Figure lla:8lock Diagram for I nstantaneous Velocity Measurements 
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Figure Ub:Mul ti-sensor Anemometer 
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' ■■■, Tangential Location-One Division=2 i t Degrees 

Figure 14c:Blade to Blade Distribution of Turbulence Intensities 
R=0. 752 ,Stat ion 2 a 
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Tangential Location-One Division=2 1 t Degrees 

Figure I5d :B!ade to Blade Distribution Of Turbulence Intensities 
T ,T .Station 2a 
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Tangential Location-One Division=2lt Degrees 

Figure l6a:Blade to Blade Distribution of Radial Vel oc i ty , R=0 . 920 
Station 2a 







w/u 
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Tangential Location-One Division=2A Degrees 

Figure l6b:Blade to Blade Distribution of Relative, Absolute Tangential 
and Axial Veloc i t ies ,R=0. 920 ,Stat ion 2a 
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Tangential Location-One Division=2 i t Degrees 

Figure l6d:Blade to Blade Distribution of Turbulence Intensities 
T ,T ,R=0. 920, Station 2a 


4 bladed inducer 3 bladed inducer 
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d - Circumferential position in degrees 

Figure 17a: Blade toBlade Variation of Stagnation and Static Heads, 

Resultant Velocity, of the Relative Flow at Station 2,R=.974 


W/U from hot-wire measurements 



0 - Circumferential position in degrees 
Figure 17b: Blade to Blade Variation of Stagnation and Static Heads, 

Resultant Velocity, of the Relative Flow at Station 2,R=.920 






k bladed -inducer • . . ' >• 3,-bladed inducer 
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Figure 17d: Blade to Blade Variation of Stagnation and Static Heads, . 

Resultant Velocity, of the Relative Flow at Station 2,R=.808 


W/U from hot-wire measurements 



Uk 


0 - Circumferential position in degrees 

Figure 17e: Blade to Blade Variation of Stagnation and Static Heads, 

Resultant Velocity.of the Relative Flow at Station 2,R=.754 



- Circumferential position in degrees 

Blade to Blade Variation of Stagnation and Static Heads, 
Resultant Velocity, of the Relative Flow at Station 2,R=.720 
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0 - Circumferential position in degrees 
Figure 17 g ; Blade to Blade Variation of Stagnation and Static Heads, 

Resultant, Velocity , of the Relative Flow at Station 2,R=.672 



W/U from hot-wire measurements 
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0 - Circumferential position in degrees 

Figure 17^ : Blade to Blade Variation of Stagnation and Stat ic 'Heads , 

Resultant Velocity, of the Relative Flow at Station 2,R=.6l7 




d Static Heads, 
t Station 2,R=.562 



A bladed inducer 3 bladed inducer 

* W/l) 
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Figure 1/ j: Blade to Blade Variation of Stagnation and Static Heads, 

Resultant Velocity, of the Relative Flow at Station 2,R=.506 
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Figure 21 : Annulus Wall Static Pressure Distribution 
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